ASYMPTOTIC EXPANSION OF THE INTEGRATED DENSITY OF 
STATES OF A TWO-DIMENSIONAL PERIODIC SCHRODINGER 

OPERATOR 



LEONID PARNOVSKI & ROMAN SHTERENBERG 

Abstract. We prove the complete asymptotic expansion of the integrated density of 
states of a two-dimensional Schrodinger operator with a smooth periodic potential. 



1. Introduction 
Let H be a Schrodinger operator 
(1.1) H=-A + V 

acting in ~R d . The potential V is assumed to be infinitely smooth and periodic with 
r cl d being its lattice of periods. We denote by = K 2 /T the fundamental quotient 
of T and by v the L°°-norm of V. We also denote by a dual lattice to T and put 
QT = R 2 /r"l\ Denote by N(X) the (integrated) density of states of the operator H. The 
density of states is defined by the formula 

Here, Hp is the restriction of H to the cube [0, L] d with the Dirichlet boundary condi- 
tions, and iV(A; A) is the counting function of the discrete spectrum of (a bounded below 
operator with compact resolvent) A. If we denote by iV (A) the density of states of the 
unperturbed operator Ho = —A, one can easily see that for positive A one has 

(1.3) iVo(A) = j^y d w d \ d/2 , 
where 

(1.4) w d 



71 



d/2 



r(l + d/2) 



is the volume of the unit ball in ~R d . There is a long-standing conjecture that for large A 
the density of states of the perturbed operator enjoys the following asymptotic behaviour 
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as A — > oo: 

oo 

(1.5) #(A)~# (A)(l + j>A-'"), 

3=1 

meaning that for each K G N one has 

(1.6) N(X) = N (\)(l + J2 e i X ~ j ) + R xM 

3=1 

with i?i<-(A) = o(A^ _K ). In these formulas, are real numbers which depend on the 
potential V. They can be calculated relatively easily using the heat kernel invariants 
(computed in [2]); they are equal to certain integrals of the potential V and its derivatives. 
Indeed, in the paper [TJ, all these coefficients were computed; in particular, it turned out 
that if d is even, then ej vanish whenever j > d/2. 

So far, formula (11.51) has been proved only in the case d — 1 in the paper [TTJ. In 
the multidimensional case, only partial results are known, see PQ, |3], [5], [9], p2], [13] . 
In particular, in [T3] it was shown that when d = 2 formula (II. 6p is valid with K = 2 
and -R(A) = 0(A~s +e ) for any positive e; in [1] it was shown that when d > 3 formula 
(11.61) is valid with K = 1 and -R(A) = 0(A~ 5 ) with some small S when d = 3 and 
i2(A) = 0(A^ In A) when d > 3. 

The aim of this paper is to establish the complete asymptotic formula (II. 5p in the 
2-dimensional case. Namely, we will prove that if d = 2, we have: 

(1.7) N(X) = -j— (A -b) + 0{\~ K ) 

47T 

for each K e N as A — > oo with 

(1.8) 6 := — ^- / ^dr. 

vol(O) y 

Note that in view of [7], it is enough to establish that (I1.6P holds for each K with some 
constants e^; then (I1.7P will follow automatically. Moreover, suppose that we have proved 
the following asymptotic formula: 

2K 2K 

(1.9) N(X) = N {X)(l +J2 e 3 X ~ j/2 + X)e i A-^ 2 ln(A)) + o(A^^). 

j=l 3=1 

Then, applying the same arguments as in [7J , together with some straightforward calcu- 
lations (one needs to compute the Laplace transform of X a In A), it is easy to show that 
(11.91) still implies (11. 7ft . Therefore, our aim will be to prove (II. 9ft . It was quite surprising 
for us when we were performing the calculations that the terms containing logarithms 
were actually 'present' in the asymptotics of N, although the coefficients Cj in front of 
these terms turned out to be zero. 
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Remark 1.1. The coefficients ij in front of logarithmic terms can be non-zero if one 
allows non-local pseudo-differential perturbations V. For example, suppose, Y = 1? and 
V is a pseudo-differential operator of order zero with the following symbol: 

(1.10) «(x,£) = [cos(27ra;i) + cos(2ttx 2 ) + cos(27t(x! - z 2 ))]Xi(|£|)X2(arg(£)), 

where xi is a smooth cut-off to the interval [l,+oo), and \i is a smooth cut-off to 
[— 0.1,7r/4]. Then the formula ( 11. 9ft is still valid, with £4 7^ 0. This can be seen by 
repeating the arguments of our paper for non-local operators and a careful computation 
of all coefficients. Since in our paper we do not consider non-local perturbations, we will 
not go into more details, but we may return to this example in a further publication. 

The method we apply to establish (11.91) consists of two parts. The first part is, es- 
sentially, the method used in [8] in order to prove the Bethe-Sommerfeld conjecture in 
all dimensions, while the second part consists of a detailed analysis of the eigenvalues 
coming from the different zones (resonance and non-resonance ones); when working in 
the resonance regions, we use some arguments from the theory of analytic functions of 
several complex variables. Dealing with the resonance regions is the part of the proof 
which at the moment we cannot extend to higher dimensions (more on this later). Now 
let us discuss the general strategy of the proof in detail (but still on a not too formal 
level) 

The first step of the proof, as usual, consists of performing the Floquet-Bloch decom- 
position to our operator (II. lft : 

(1.11) H = [ H(k)dk, 

J® 

where H(k) = Hq + V(x) is the family of 'twisted' operators with the same symbol as 
H acting in 

f):=L 2 (0). 

These auxiliary operators are labelled by the quasi- momentum k G 0^; the domain D(k) 
of H(k) consists of functions / G H 2 (0) which are restrictions of functions / G Hf oc (M. 2 ) 
satisfying the following condition: /("y + x) = e* k ' 7 /(x), 7 G T. We refer the reader to 
[TO"] for more details about this decomposition. Now it would be useful to introduce a 
different density of states 

(1.12) N(X) := [ N(X,H(k))dk 

Jot 

which is more convenient to deal with. It is known (see e.g. [TU]) that 

(1.13) N(X) = ^iV(A). 

Therefore, for our purposes it would be enough to prove (jl.9p for iV instead of N. 

Note that we can assume without any loss of generality that J V(x)<ix = 0. Indeed, 
otherwise we consider a new operator Hi := H — b (b is defined in ( II .81) ). Note that 
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Hi = —A + V\ and the constant Fourier coefficient of V\ := V — b vanishes. Since 
N(X;H) = N(X — b,Hi), we see that asymptotic formulas (11.71) for H and Hi are 
equivalent. Therefore, we can (and will) always assume that / V(x)dx = 0. 

Next, instead of trying to prove (11.91) for all values of A, we will prove it assuming 
that A is inside a fixed interval: A G [A n , 16A„], where A n = 4 n A is a large number, and 
we will allow the coefficients in (11.91) to depend on n, although the remainder should be 
uniform in n. In Section 3, we will show that if we can prove these asymptotic formulae 
for all n with coefficients growing not too fast, this would imply the validity (II. 9p for all 
A. The reason we require this reduction is the following: on later stages of the proof, we 
will decompose the phase space (i.e. the space where the dual variable £ lives) into two 
regions: resonant and non-resonant zones. The resonant zones are, roughly speaking, 
the strips of some width a. The value of a cannot be chosen the same for all values of A, 
since we need a to be of order A 1//6 . Thus, when we increase A, at some stage we will have 
to increase the value of a, and this can result in changes of the asymptotic coefficients 
in (II. 9p . However, if A runs over a fixed interval [A n , 16A n ], we can keep a fixed and thus 
the coefficients of our asymptotic expansion (11.91) will be fixed as well. Thus, starting 
from section 4, we will be assuming that A G [A n , 16A n ] and n is fixed. 

The next step is to assume that the potential V is a finite trigonometric polynomial 
whose Fourier coefficients 

(1.14) V (m) := -=L= [ ^(x)e~^x, m G I* 

a/vo1(0) Jo 

vanish when |m| > R. More precisely, we replace the original potential 

(1.15) V(x) = -7^== £ V{m)e— 

v vol (°) S 

by the truncated potential 

(1.16) V(x)= J— n™)e- imx , 

V vol (0) meB(i?)nr t 

where B(R) is a ball of radius R centered at the origin. Here, R = R n is a parameter 
which grows as a small positive power of A n (for example, R n = A^ 48 ). It is easy to 
justify he fact that the error introduced by changing the potential in such a way is small; 
this is where we use the fact that the original potential is infinitely smooth. However, 
this truncation leaves us with an additional tedious job of checking how all the important 
estimates depend on R n . 

Next, our aim is to construct a good approximation of all the eigenvalues of all oper- 
ators H(k) simultaneously (to be precise, we will need to approximate only eigenvalues 
which are inside the interval [A n — 100f , 16A n + 100w]). In Section 4, we discuss what 
exactly we mean by such a simultaneous approximation and prove that, if this approxima- 
tion satisfies a bunch of additional properties (in particular, the approximating function 
needs to behave in a proper way in a specially chosen coordinate system which should 
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also satisfy certain properties), then asymptotic formula ( 11. 9ft would follow automati- 
cally. This is done in Lemma 14.101 Unfortunately, we will not be able to use this lemma 
without modifications further on in the paper, but at least this lemma (and the proof of 
it) shows us which properties we are aiming for. 

The main part of the paper, Sections 5-7, is devoted to the construction of an approxi- 
mation of all the eigenvalues of all operators H(k); the existence of such an approximation 
was assumed in Section 4. The main tool during this construction will be an abstract 
result from perturbation theory - Lemma 15.11 This lemma allows us, under certain 
conditions, to study the spectrum of an operator PHP instead of the spectrum of an 
operator H. Here, H = H + V is a bounded below operator with compact resolvent, V 
is bounded, and P is a spectral projection of H . Since the formulation of this lemma is 
rather involved, let us illustrate what it says by considering a special case. Assume that 
P is a further sum of spectral projections of H , P = J2j=o Pj suc h that the matrix of V 
in the basis corresponding to Pq, Pi, . . . , Pj, Pj+i := / — P is block-three-diagonal (i.e. 
PjVPt = whenever \j — 1\ > 1). Assume also that A = X(H) is an eigenvalue of H and 
that the distance from the spectra of PjH Pj (j = 1, . . . , J+ 1) to A is at least a, where a 
is sufficiently large, so that P HP is 'essentially responsible' for the eigenvalue A. Then 
the operator PHP has an eigenvalue A' such that |A — A'| <C a~ 2J . In applications, a 
will be of order A 1//g , so by choosing sufficiently large J, we can make our approximation 
as precise as we wish. 

We are going to apply Lemma 15.11 by constructing various projections P such that 
the operator PHP has an eigenvalue close to an eigenvalue of H. Roughly speaking, 
each point £ from the phase space such that |£| 2 is close to A generates such a pro- 
jection P = P(£). The structure of P($) depends on the exact location of £ in the 
phase space. There are two types of points resonant and non-resonant ones. For 
non-resonant points the structure of P(£) is relatively simple, the operator PHP has 
a unique eigenvalue close to A, and we can find this eigenvalue using the standard ap- 
proximating procedure (for example, the Banach contraction mapping theorem). Having 
constructed this approximation, we are ready to start computing N(X); it is a relatively 
straightforward (but slightly tedious) task to compute the contribution to the density 
of states coming from the non-resonant regions. The logarithmic terms appear on this 
stage (resonant regions do not produce any logarithms). 

In the case of the resonant the structure of P(£) is more complicated, and therefore 
it is much more difficult to compute a contribution to the density of states coming from 
the resonance zones. The main problem lies in the fact that the approximation formula 
for the resonant eigenvalues is not explicit: it expresses eigenvalues of PHP in terms of 
the eigenvalues of an expression A + eB, where A and B are explicitly given symmetric 
matrices and e ~ is a small parameter which also depends on £ in an explicit way. 
Of course, one can expand the eigenvalues of A + eB in powers of e, but the coefficients 
in this expansion will not be uniformly bounded in £, so we will not be able to integrate 
this expansion in Thus, we need to analyse the situation deeper. Let us denote 
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by P the projection onto the kernel of A. (We are interested in the perturbation of 
zero eigenvalues of A.) Then a priori there are two reasons why the coefficients in the 
asymptotic expansion of A + eB can be large: either A has eigenvalues close to zero 
(not the case in our situation), or the operator PBP has eigenvalues close to each other. 
The latter possibility is actually occurring in our problem. However, it turns out that 
PBP is 'essentially' unitary equivalent to a one-dimensional Schrodinger operator with 
quasi-periodic boundary conditions on the interval. Therefore, there could be no more 
than two eigenvalues of this operator located near each other (at this place we strongly 
use the fact that our operator H is two-dimensional). The rest of the computations is 
similar to the non-resonance regions, only instead of solving equation /x + G(/i) = A 1 / 2 
like we did in the non-resonance region (and where we used implicit function theorem), 
now we have to solve the equation /i 2 + Xi(/i)/i + X 2 (/i) = 0. The tool for dealing with 
equations of this type comes from the theory of functions of several complex variables 
and is called the Weierstrass Preparation Theorem. After using this theorem, we obtain 
the expressions for eigenvalues in the non-resonance regions; these expressions are no 
longer analytic in A 1//2 , but contain square roots of analytic functions; however, these 
square roots will cancel after integration in £ to produce an asymptotic formula which 
contains only powers of A 1 / 2 . 

The rest of the paper is organised as follows: in the next section we give all necessary 
definitions and basic facts (the Weierstrass Preparation Theorem and corollaries from 
it). In Section 3 we reduce the problem of finding an asymptotic formula valid for all A to 
the problem of finding such a formula valid only for A inside a fixed interval. In Section 
4 we describe what exactly we mean by a simultaneous approximation of all eigenvalues 
of all H(k) and give some idea about the general strategy of the proof. In Section 5 
we formulate auxiliary results which were proved in [5] and introduce the partition of 
the £-plane into resonance and non-resonance regions. In Section 6 we deal with the 
non-resonance regions, and, finally, in Section 7 (the most complicated one) we compute 
the contribution to the density of states from the resonance zones. 
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author was partially supported by the LMS grant. Both of us are very grateful to Yu. 
Karpeshina and A. Sobolev for useful discussions. Finally, we thank the referees for 
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2. Notation and basic facts 

Let T be a lattice in M 2 . We denote by = M 2 /r the fundamental domain of T, by 
r^" the lattice dual to T, and by 0^ = M 2 /rt its fundamental domain. 
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For each vector x e M 2 we denote by x -1 the result of rotation of x by — | and 
n(x) := a, assuming x / 0. If Xi,x 2 G M 2 are two non-zero vectors, we denote by 

0(xi,X2) the angle between them (0 < <fi < it). 

If ^ G M 2 , there exists unique decomposition £ = 7 + k with 7 G and k 6 We 
call 7 =: [£] and k =: {£} resp. the integer part and the fractional part of 

If is a bounded below self-adjoint operator with compact resolvent, then Hj{H) is 
its j-th eigenvalue (counting multiplicities). 

We also assume that the average of V over is zero. 

By C or c we denote positive constants, The exact value of which can be different each 
time they occur in the text, possibly even each time they occur in the same formula. On 
the other hand, the constants which are labeled (like Ci, C3, etc) have their values being 
fixed throughout the text. Given two positive functions / and g, we say that / ^> g, or 
g <C /, or g — 0(f) if the ratio j is bounded. We say / x g if / g and / <C g. 

The results in the rest of this section are quoted from [3]. 

Theorem 2.1. (The Weierstrass preparation theorem). Let F be analytic and bounded 
in a neighborhood uofOin C n and assume that F(0, z n )/ 'z p is analytic and ^OatO. (In 
other words it means that (d^F/dz^)(Q, 0) = for j = 1, . . . ,p— 1, and (d p F/dzP)(0, 0) ^ 
0). Then one can find a polydisc D C w such that every G which is analytic and bounded 
in D can be written in the form 

(2.1) G = qF + r, 

where q and r are analytic in ~D, r is a polynomial in z n of degree < p (with coefficients 
depending on z' = (z±, . . . , z n _i) ) and 

(2.2) sup |g| < Csup|G|. 

D D 

The representation is unique. 

Remark 2.2. It follows from the proof that polydisc D and constant C can be chosen 
to depend only on sup \F\, p and ((d p F/dz^)(0, 0)) _1 . 

Now, choose G := z p and put W := z p — r, h := q^ 1 . We have 

Corollary 2.3. If F satisfies the hypothesis of Theorem \2.1\ then one can write F in a 
unique way in the form 

(2.3) F = hW, 

where h and W are analytic in a neighborhood uj 1 ofO, h(0) 7^ ; and W is a Weierstrass 
polynomial, that is, 

p-i 

(2.4) W(z) = < + 5>(/)4, 

3=0 

where aj are analytic functions in a neighborhood of vanishing when z' = 0. 
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Moreover, 

(2.5) supfl/il + l/il-^+suplWn <Ci, 

and to' and C\ depend only on sup \F\, p and ((d p F/dz^)(0, 0)) _1 . 

Corollary 2.4. Assume that a set of functions F := {F} satisfies the following proper- 
ties: 

1. Functions F are analytic in a neighborhood u of 0. 

2. For somep we have (d j F/dzl)(0, 0) = for j = 1, . . . ,p-l, and (<9 p F/<9^)(0, 0) ^ 

0. 

3. We have the bounds 

(2.6) sup sup |F| < C, sup 1(^/30(0,0) | - 1 < c - 

Then there exist a neighborhood to' and a constant C\ such that for any F e F the 
representation (I2.3P and estimate (12.51) ZioZd. Moreover, u' and C\ are uniform with 
respect to {F} and depend only on u, p and constant C from (12.61) . 

3. Reduction to a finite interval of spectral parameter 

The main result of our paper is the following theorem (or, rather, the corollary from 
it; we put p := vA): 

Theorem 3.1. For each K G N we have: 

K K 

(3.1) N(p 2 ) = np 2 + J2 ejP~ j + l*p£%^ + o(p 

3=0 j=2 

as p — ► oo. 

Once the theorem is proved, it immediately implies 
Corollary 3.2. For each K G N we have: 

(3.2) N(X) = -L A - -Lj ^ ^(x)rfx + 0(\-K) 
as A — > oo. 

Proof. First of all, we notice that [2] implies that 



-K\ 



/"OO °° 

(3.3) / e- tx N(X)d\ ~ r^/ 2 Qjt j 

as t — > 0+, where q 3 - are constants depending on the potential. Now the corollary follows 
from theorem 13.11 property fl!.13p . and calculations similar to that of [7j. Indeed, [7] 
implies that if all coefficients ej vanish, then all coefficients ej, j > 0, vanish as well. 
It remains to show that all coefficients ej vanish. Suppose, this is not the case. We 
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consider separately even and odd values of j. Suppose first that e^k is the first non-zero 
even coefficient with hats. Then we consider the following integral: 



(3.4) I(t) := I e- xt \- k \n\d\ 



and, after elementary calculations, find that the asymptotic expansion of I(t) as t —>■ 0+ 
contains a term In 2 t with a non-zero coefficient. This term is absent in the Laplace 
transform of other terms from the expansion (13. f I) . Thus, our assumption that e^k 7^ 
contradicts (13.31) . 

Suppose now that &2k+i is the first non-zero odd coefficient with hats. Then, similarly 
to the previous case, we consider the following integral: 



(3.5) I(t) := / e- A 'A-( 2fc+1 >/ 2 lnAciA 



oo 



and find that the asymptotic expansion of I(t) as t — > 0+ contains a term t( 2fc_1 )/ 2 ln£ 
with a non-zero coefficient. This term is absent from the Laplace transform of other 
terms from the expansion (13. f p . Once again, we have reached a contradiction with (13.31) . 
Thus, all coefficients e m vanish, and our corollary follows from [7]. 

□ 

The rest of the paper is devoted to proving Theorem 13. f I 

To begin with, we choose sufficiently large p > 1 (to be fixed later on) and put 
p n = 2p n _x = 2 n p ; we also define the interval = [p„,4p n ]. The proof of the main 
theorem will be based on the following lemma: 

Lemma 3.3. For each M e N and p G I n we have: 

6M 6M 

(3.6) N(p) = np 2 + Y,e J (n)p-i + lnp ^(n)p^ + 0(p~ M ). 

3=0 j=2 

Here, ej(n), ij(n) are some real numbers depending on j and n (and M) satisfying 

(3.7) ej{n) = O(p^), e» = O(pf^). 

The constants in the O -terms do not depend on n (but they may depend on M). 

Remark 3.4. Note that (13.61) is not a 'proper' asymptotic formula, since the coefficients 
ej(n) are allowed to grow with n (and, therefore, with p). 

Let us prove theorem 13.11 assuming that we have proved lemma 13.31 Let M be fixed. 
Denote 

6M 6M 

(3.8) N n (p 2 ) :=7rp 2 + ^e>)p- J +lnp^e>)p^. 
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Then whenever p £ J n := J n _i D I n = [p n , 2p n ], we have: 

6M 6M 

(3.9) N n (p 2 ) - N n ^(p 2 ) = 5>(n)p^' + hip ^(n)p-^ 
where 

(3.10) ^( ri ) := e j( n ) ~ e j( n ~ 1)) ij( n ) := — Cj(n — 1). 

On the other hand, since for p £ J n we have both iV(p) = N n (p) + 0(p~ M ) and iV(p) = 
^n-i(p) + 0(p- M ), this implies that W/ 3 " 3 ' + ln P E;f 2 *i(™)P~ j = 0(p~ M ). 

Claim 3.5. For eac/i j = 0,...,6M we /lave: tj(n) = 0(pP~ M In p n ) and £j(n) = 

o(pr M )- 

Proof. Put a; := p _1 . Then Ej=o ^j( n ) x ^ — hix E^=2 ^:/( n ) a> ' = 0(Pn M ) whenever x E 
Then 

6M 6M 

(3.11) P(y) := 5>iW + E^W ln ^ = °W 

whenever y e [|, 1]. Consider the following 12M functions: y j (j = 0, ...,6M) and 
y^lny (j = 2,...,6M) and label them h\(y), ...h\2M{y)- These functions are linearly- 
independent on the interval [|,1]. Therefore, there exist points yi,...,yi2M £ [§, 1] 
such that the determinant of the matrix (hj(yi))^ 2 ^ l is non-zero. Now (13.111) and the 
Cramer's Rule imply that for each j the values Tj(n) and fj(n) are fractions with a 
bounded expression in the numerator and a fixed non-zero number in the denominator. 
Therefore, Tj(n) = 0(1) and fj(n) = 0(1). This shows first that tj(n) = 0(p^ M ) and 
then that tj(n) = 0(p!~ M In p n ) as claimed. □ 

Thus, for j < M, the series Em=oA?'( m ) * s absolutely convergent; moreover, for such 
j we have: 
(3.12) 



n 



e 3 [n) 



e;(0) + $>H = e,(0) + ^ t,(m) + 0(pT M In p n ) =: Cj + 0(pT M In p n ), 

m=l m=l 

where we have denoted Cj := e,-(0) + Em=i tj( m )- Similarly, for j < M we have 

n oo 

(3.13) e» = e,(0) + £ t,(m) = e,(0) + £ t,(m) + 0(p^ M ) =: e 3 - + 0(fP- M ), 

m=l m=l 

where we have denoted ij := ij(0) + Em=i tj( m )- 
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4j+7 

Since ej(n) = 0(p n 5 ) (it was one of the assumptions of lemma), we have: 

6M 

7_M M 

(3-14) \ e o( n M j = 0(pV 5 ) = 0( Pn 6 ), 

j=M 

assuming as we can without loss of generality that M is sufficiently large. The sum with 
hats on is estimated similarly. Thus, when p G I n , we have: 

M-l M-1 

(3.15) N(p) = up 2 + e iP~ J + ^P~ j ln P + 0(p- M \np) + O(p-f). 

3=0 3=2 

Since constants in O do not depend on n, for all p > p we have: 

M-1 M-1 

^ ' ^ [AT/6] [A//6] 

= np 2 + ^ e j p~- J + ^2 ejp~ J lnp + 0(p"~). 

j=0 j=2 

Taking M = + 1, we obtain fl3TT|) . 

The rest of the paper is devoted to proving lemma 13.31 

4. Description of the approach. Integration in new coordinates 

From the previous section it is clear that we can study the density of states N[p) as- 
suming that p G I n . Throughout the paper we will assume that n is fixed and sometimes 
will omit index n from the notation; however, we will carefully follow how all estimates 
depend on n. If we need to make sure that p n is sufficiently large, we will achieve this 
by increasing p , keeping n fixed. 

First, we discuss the general strategy. In this section we describe how to construct the 
asymptotic formula for N(p) using certain objects (mappings / and g and coordinates 
(r, $) satisfying certain properties); in the next sections, we will construct these objects. 

Let us fix sufficiently large n, X = p 2 with p 6 7 n , and denote 

(4.1) A = A(p) := {$, E K 2 , HI 2 E [A - 100v, A + lOOu]}, 

where v := ||^||oo- Obviously, A is an annulus of width ~ p _1 . We also fix a number 
MgN. Our aim is to construct good approximation of the eigenvalues lying close to A. 
Namely, we will construct two mappings /, g : M 2 — > R such that for each £, /(£) is an 
eigenvalue of If ({£}); moreover, / : {£ G IR 2 , {£} = k} — > a(H(k)) is a bijection for each 
k (here, we count all eigenvalues of H(k) according to their multiplicities; the functions 
f,g depend on n, M and p). The difference — g(£)\ is required to be sufficiently 

small at least when £ G A, namely, we postulate that the following two properties hold: 

« \f(Z)-9m<Pn M toT$eA; 

(h) 1/(1) - l£ 2 || < 2v, similarly, |^(^) - |^ 2 || < 2v. 
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Notice that the second property implies that if £ G" A, then the following three in- 
equalities are equivalent: /(£) < A if and only if g(£) < A, and this in turn happens if 
and only if |£| < p. 

Remark 4.1. Rigorously speaking, the functions we will construct will satisfy property 
(i) not in the whole annulus A, but in a slightly smaller annulus G R 2 , |£| 2 G [A — 
90v , A + 90f]}. Indeed, in the process of constructing / and g we will have to reduce 
the width of the set A by 2v several times. One obvious solution to this problem would 
be to introduce sets Ai = {$, G R 2 , |£| 2 G [A - 98v , A + 98t>]}, A 2 , etc. However, this 
would introduce extra notational complexity to a paper which is already overburdened 
with notation. Thus, we will keep calling A all annuli of slightly smaller width whenever 
necessary. 

Finally, we will construct function g in such a way that it satisfies some asymptotic 
formulas. The next lemmas describe why these functions are going to be useful. Denote 

B f (X) := r\(-oo,X]). 

Lemma 4.2. Suppose f : R 2 —>■ R is a measurable mapping such that f : {£ G 
R 2 , {£} = k} —>■ cr(H(k)) is a bijection (including multiplicities) for each k. Then 
N(X)=v6L(B f (X)). 

Proof. Denote by XB f (\) the characteristic function of Bf(X). By Fubini's theorem we 
have: 



vol(5 / (A))= / XB fi x)(m 
(4.2) Jw 



#{£,{£} = k&/(£)<A}dk= / iV(A,k)rfk = iV(A). 
ot Jot 



□ 



Our next task is two-fold: to show that under certain conditions we can replace Bf in 
lemma H~2l by B g so that the error is not too big and, secondly, to compute vol(B g (X)) 
(or, at least, to expand this volume in powers of A). Unfortunately, assumptions (i) 
and (ii) on functions / and g made above are not the only necessary requirements to do 
this job: we also need to check that function g behaves in a 'nice' way in some suitable 
coordinates. Since the complete set of required conditions looks rather nasty, we will 
introduce these conditions slowly, on at a time, to show why each particular condition is 
required. First, we check that the polar coordinates could do the trick. 

Lemma 4.3. Let A = p 2 be fixed and let N be a fixed natural number. Suppose f : R 2 — ► 
R is a measurable mapping such that f : G R 2 , {£} = k} — > a(H(k)) is a bijection 
for each k (counting multiplicities). Suppose, g : R 2 — > R is a measurable mapping and 
that f,g satisfy properties (i) and (ii) above. Suppose also that ^(re 4 *) ^> p whenever 



re 



** G A. Then N(X) = vol(B g (X)) + 0(p~ A/ ) 
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Proof. Assumptions of lemma (namely, property (ii) above) imply that the symmetric 
difference Bf(X)AB g (X) C A. The boundary of the 'ball' B g (X) is a subset of A; since 
the function g = g(r,<p) is increasing in r, for any fixed 4>q the intersection of B g (X) 
with any semi-infinite interval {re 1 * , r £ [0, oo)} is an interval {re 1 * , r £ [0, Z]}, where 
Z = Z((j)o) is a well-defined function. Since |^ > Cip, we also have that if r < Z(<f>) — 
G^ x p- M ' x , then X-g(re i(f> ) > p~ M , and so /(re**) < X and £ = re** £ £/(A). Similarly, 
if r > Z((j)) + C^p- 1 *- 1 , then #(re^) - A > p~ M , and so /(re**) > A and $, = re** & 
B f (X). Thus, the symmetric difference B f (X)AB g (X) C {£ = re 1 * £ M 2 , \r - Z(<f>)\ < 
Cx 1 p~ M ~ 1 } and thus vol(S/(A)AS ff (A)) < p~ M . Together with lemmalOl this finishes 
the proof. □ 

Later on, we will apply lemma 14.31 in a more general situation, when r is not precisely 
the radial coordinate, but 'close' to the radial coordinate in a certain sense; more pre- 
cisely, we will need the following statement (with proof being exactly the same as proof 
of lemma I 



Corollary 4.4. Let § be a curve of length <C I, and let (r, $) (r £ 1R + , $ £ §) be 
coordinates in A(p) such that the Jacobian | §|^y |<S p- Suppose, f : M 2 — > K is a 
measurable mapping such that f : £ M 2 , {£} = k} — > o~(H(k)) is a bisection for 
each k (counting multiplicities). Suppose, g : M? — > R a measurable mapping and 
that f,g satisfy properties (i) and (ii) above. Suppose also that |^(r, $) ^> p. T/ien 
iV(A)=vol( J B 9 (A)) + 0(p- M ). 

Remark 4.5. Obviously, lemma H~3l is a special case of corollary 14.41 with S = S 1 being 
a circle of radius 1 centered at the origin and (r, $) being the usual polar coordinates. 

Remark 4.6. Suppose that another set of coordinates (f, $) satisfy slightly different 
conditions: $ £ S, where S is a curse of length <C p, but the Jacobian I fff'ju |<C 1. 



Then the coordinates (r, $) := (f, — ) satisfy all assumptions of corollary I4.4[ so the 
conclusion of this corollary will also be valid for such coordinates. We will be using both 
types of coordinates, depending upon convenience. 

Remark 4.7. The coordinates (r, <£>) which we will introduce in further sections will be 
defined simultaneously for all p £ J n , i.e. they will be defined for all points 

(4.3) ££A^ :=U pGln A(p). 

Unfortunately, we need to make our assumptions about the coordinate system (r, $) 
even more complicated. First of all, we will need to use different coordinates systems in 
different parts of A^ n \ so we assume that we have a decomposition of A^ as a disjoint 
union: 

(4.4) A {n) = uf =1 .A[ n) ; 

for simplicity, we assume that all sets Af 1 are open, and treat (I4.4p modulo points on 
the boundaries of these sets. We also assume that there is a coordinate system (r, $) 
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in each (r(£) G M + , $(£) G §/, where §; is a curve of length <C 1) and that this 
system satisfies all assumptions of corollary 14.41 Whenever we talk about the Jacobian 
| |, we will assume that it is defined only at points £ located inside some A\ , 

i.e. the Jacobian is not defined for points on the boundary of A\ n \ Other conditions we 
always assume are: r(£) ~ |£|, and for each fixed $o the intersection 

(4.5) A[ n) n{£ = (r, $ ), r G [0, oo)} = = (r, $ ), r G [n, r 2 ]} 

is an interval with endpoints £ x = (r 1; $ ) an d £2 = ( r 2; ^0) satisfying |^ x | 2 = p 2 n — 100f 
and I £ 2 1 2 — (4p n ) 2 + 100u. The latter condition, while looking rather horrific, is easy to 
check and will be always automatically satisfied in our constructions. Roughly speaking, 
it is needed to ensure that the curve {£ = (r, <3> ), r G [0, 00)} (which happens to be a 
semi- infinite interval in all our constructions) cannot enter or leave A ™ from the 'sides'. 
Technically, it is required to make sure that formulas ( 14. 18ft and ( 14. 19ft imply ( 14.201) . 

Let us introduce more notation. Put 

(4.6) A + := {Z G R 2 , g{t) < p 2 < \Z\ 2 } 
and 

(4.7) A- := {£ G R 2 , |£| 2 < p 2 < g(£)}. 
Lemma 4.8. 

(4.8) vol(B g (p 2 )) = up 2 + voli+ - voli - . 

Proof. We obviously have B g (p 2 ) = B(p 2 )uA + \A- . Since Ar C B(p 2 ) and A + C]B{p 2 ) = 
0, this implies (Ojl . □ 

Remark 4.9. Property (ii) of the mapping g implies that we have A + , A~ C .A. Thus, 
statements I4.2H4.8I imply that in order to compute N(X), we need to analyse the be- 
haviour of g only inside A. 

In order to apply corollary 14.41 and lemma H~8l for computing the asymptotic behaviour 
of N(p), we need even more assumptions. Roughly speaking, the next lemma says that 
whenever all objects involved enjoy a power asymptotics at infinity, then so does B g (X). 

Lemma 4.10. Let I G N and a G (0, 1) be fixed. Suppose that all assumptions of 

corollary \4.4\ and remark \4^\ are satisfied and that for fixed <3> the point £ = (r, $) G A^ 
has an absolute value |£| which has an asymptotic expansion in powers of r: 



(4.9) HI = r(£) 



J 



+ 0(r(Z)- M ) 
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d\t\ 
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(4.10) 



Or 



i=2 



9(t) = r(& 



Suppose also that the function g enjoys the following asymptotic behaviour in r(£) when 

\ 

+ 0(r($)- M ). 

Finally, suppose that the Jacobian also satisfies an asymptotic formula: 
d(x,y) 



(4.11) 



( |" M+2 j 



(4.12) 



d(r, $) 



-Af\ 



i=i 



All functions aj, dj, etc. are measurable and bounded (but not necessarily continuous) 
functions of $ and are 0(p aj ), a < 1. Taen 



(4.13) 



voi(i + n yi}" ) - voi(i- n Af n) ) = p 2 E 6 ^~ j + °(p~ m ) 

3=1 



and all bj are 0(p aj ). 

Remark 4.11. It may seem strange that absolute value of the power in the remainder 
term in the above formulas is smaller than the upper summation limit. This is caused by 
the fact that the coefficients aj, dj, bj, etc. are allowed to grow together with p: compare 
this with remark 13.41 

Proof. First of all we notice that without loss of generality we can assume that 

([M±2] \ 
i + E tiWtMt)-* 
J 

since corollary 14.41 implies that the error caused by using this approximation is 0(p~ M ). 
Let us for a moment fix some value $o- Then the RHS of f)4.14p is an increasing function 
of r for sufficiently large r. Let us call by Q 1 = Q\ the inverse function to f)4.14p . i.e. 

\ 



(4.15) 



{Q\t)f 



( [t=£\ 

\ j=1 



t. 
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It is an easy exercise to show that the function Q 1 also enjoys the asymptotic behaviour 
as |£| — > oo: 

(4.16) 



Ql(t) = t 



1/2 



/ [M±l] \ 
3=1 



0(t 



-M/2\ 



and that the coefficients bj = 0{p a ^). Note that Q l is also monotone increasing, so the 
inequality g(£) < p 2 is equivalent to r(£) < Q\(p 2 ). 

Equation (14.101) implies that the RHS of (14. 9 p is an increasing function of r. Let us 
denote by Q 2 = Q\ the inverse function to it. Then again it is easy to show that Q 2 also 
enjoys the asymptotic behaviour: 



(4.17) 



Q\t)=t 



\ 



i + J2 hm))t 

3=1 



-3 



0(t 



J 



with bj = 0(p aj ). Moreover, Q 2 (t) is a monotone function for large £, so the inequality 
|^| < p is equivalent to r(£) < Q%(p). 
Now we can re- write definitions fl4.6p - fl4.7l) in the following way: 



(4.18) 
and 
(4.19) 
Therefore, 

(4.20) 



A" 



A~ 



^ 2 iQ 2 M<r(ti)<QUp 2 )} 
l 2 ,QUp 2 )<r(H)<Q 2 M}. 



voi(i+nyi[ n) )-voi(i-nyi 



(n)v 

i , 




d{x,y) 
d(r, $) 



drd§. 



Now (jUISD follows from (|Q5j) . fHTTj) . and (|Q2|) . 



□ 



Remark 4.12. When applying lemma H.10I later, we will first establish asymptotic for- 
mula (14.111) only for £ e A(p) with a fixed p G /„. After this formula is established for 
each p E I n , we just check that the coefficients do not depend on the particular choice 
of p, so this formula holds for all £ G A 



(n) 



Remark 4.13. Note that logarithms have made a brief appearance in the RHS of (I4.20p 
before being canceled out. 

Remark 4.14. Lemma [4.101 gives us only a priori estimates on coefficients bj. In fact, 
we will be able to say more about them. For example, since A + , A~ C A and vol .A 1, 
this implies that the LHS of (14.131) is bounded and, thus, leads to additional restrictions 
on the first several coefficients bj. Later we will come back to this discussion. 
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5. Abstract perturbation results and decomposition into invariant 

subspaces 

In this section, we begin the construction of the mappings /, g with properties (i), (ii) 
stated in the previous section. 

First, we formulate the abstract result which was proved in [8] (lemma 3.2 and corollary 
3.3); see introduction for an informal discussion of this result. 

Lemma 5.1. Let H$, and V be self-adjoint operators such that Hq is bounded below 
and has compact resolvent and V is bounded. Let {P 1 } (I = 0, . . . , L) be a collection of 
orthogonal projections commuting with H such that if I ^ n then P l P n = P l VP n = 0. 
Denote Q := I — ^2 P l . Suppose that each P l is a further sum of orthogonal projections 
commuting with H : P l = EjU p j such that p }V p l = for \j - 1\ > 1 and P^VQ = 
if j < ji. Let v := ||V|| and let us fix an interval J = [Ai, A 2 ] on the spectral axis which 
satisfies the following properties: spectra of the operators QH Q and PjH Pj, j > 1 lie 
outside J ; moreover, the distance from the spectrum of QH Q to J is greater than 4t> 
and the distance from the spectrum of PjH Pj (j > 1) to 3, which we denote by aL, is 
greater than 12v. Denote by fi p < ■ ■ ■ < \i q all eigenvalues of H = H + V which are 
inside J. Then the corresponding eigenvalues jl p , . . . ,jl q of the operator 

H :=J2 plHpl + Q H oQ 
i 

are eigenvalues of Y2i P l HP l , and they satisfy 

\u r — uJ < max 
i 

all other eigenvalues of H are outside the interval [Ai +v, A2 — v]. More precisely, there 
exists an injection G defined on the set of eigenvalues of the operator P l HP l (all 
eigenvalues are counted according to their multiplicities) and mapping them to the subset 
of the set of eigenvalues of H (again considered counting multiplicities) such that: 

(a) all eigenvalues of H inside J have a pre-image, 

(b) If fi r G [Ai + 2v , A2 — 2v] is an eigenvalue of Yli P l HP l , then 

\G(jjL r ) — Hr\ — m ; ax 

and 

(c) G(fi r ) = fi r+ T{H), where T is the number of eigenvalues of QHqQ which are 
smaller than X\. 

Finally, we have: \\H — H\ \ < 2v. 



.7=1 



31 
.3=1 
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Let us fix n and M, and let A = p 2 be a real number with p G J n - Consider the 
truncated potential 

(5.1) V'(x)= £ t>(m)e ra (x), 

mGS(i? n )nrt 

where 

e m (x) := -J—e'K*), mert, 
Vvol(O) 

and 

(5.2) V^m) = / y(x)e_ m (x)rfx 



Jo 

are the Fourier coefficients of V. R n is a large parameter the precise value of which will 
be chosen later; at the moment we just state that R n ~ fP n with p > being small. 
Throughout the text, we will prove various statements which will hold under conditions 
of the type R n < p P n . After each statement of this type, we will always assume, without 
possibly specifically mentioning, that these conditions are always satisfied in what follows; 
at the end, we will choose p = mm pj. 
Since V is smooth, for each m we have 

(5.3) sup|\/(x)-\/ / (x)|«^ m . 

This implies that if we denote H'(k) := H Q (k) + V with the domain 2)(k), the following 
estimate holds for all n: 

(5.4) MHQt)) - ^-(#'(k))| « Rr n m « p;r- 

Thus, if we choose sufficiently large m, namely m > M/p, we can safely work with the 
truncated operator H' instead of the original operator H . 
For each natural j we denote 

(5.5) e i := rt n B(jR n ), e := {o}, ej : = e, \ {o}. 

We also choose a number M := 3M. Each vector from <y G O' -, generates a one- 
dimensional linear space {t"y, t e R}. The intersection {i-y, t e M} D contains two 
vectors with the smallest length. We call such vectors the primitive vectors. Note that 
if 6 is a primitive vector, then so is —6. Let &i, . . . , &l be the set of all the primitive 
elements of ©g^- We choose the labeling in such a way that if we take n(0i) and start 
rotating it counterclockwise, we meet n(#2), n (^3), etc. in consecutive order. 

Lemma 5.2. If 7, u G @' 15A ^ ore two linearly independent vectors, then the angle 
0(7, v) > R~ 2 for large R n . 

Proof. It is a simple geometry (and was proved, e.g. in [S], lemma 4.2 and corollary 
4.3). ' □ 
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Corollary 5.3. Under assumptions of lemma EH we have |(n(7), n(i^- L ))| ^> R n 2 for 
large R n . 

Let = 6i be a primitive vector which we consider fixed for the moment. Let us 
introduce cartesian coordinates on a plane where the first axis goes along d ± , and the 
second axis goes along 6. We call this set of coordinates coordinates generated by 6. 
Sometimes, we will also need the cartesian coordinates which are fixed and independent 
of the choice of 6f, we will call such set of coordinates universal coordinates. 

This choice of coordinates generated by means that each £ £ M 2 has coordinates 
(^1,^2), where £1 = (£,n(0 ± )) and £2 = (£, n(0)). Let us fix this coordinate system for 
now. We also define a = a n to be the smallest real number which satisfies two conditions: 

(5.6) a>pV3 & ^_^ N . 

1/3 1/3 1/3 

In particular, we have P T <a< p 1 rI° + \e\/2 + l < 2pT. Now we can make the following 
definitions: 

(5.7) A(0):={£€M 2 , \(£,n(0))\ < a}, 

(5.8) E 1 (0) := {£ £ A(p) PI A(0), (£, ± ) > 0}. 

Obviously, the intersection PI A(0) consists of two connected components, and the 
condition (£, ) > chooses one of them. We also define 

(5.9) S 2 (0) := {rj = £ + W, £ £ Hx(0), * £ E}, 

(5.10) E 3 (0) :=H 2 (0)nA(0), 
and 

(5.11) E 4 (0):=(yt(p)nH 2 (0))\H 3 (0). 




Figure 1 
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Figure 3 

Lemma 5.4. Suppose R n < pi/ 10 , £ g A(0) andj = W G Q'^. Then ||£ + 7| 2 - |£| 2 | > 

1/3 
Pn ■ 

Proof. Indeed, we have 

ii^+7r-i^ri>i(t7)i-i7i 2 »a>py 3 . 

□ 

Lemma 5.5. Let rj G Then \rfi — p\ <C p^ 1 ^ 3 . 

Proof. Indeed, since rj G A(p), we have |r/| 2 = rjf + rfe = p 2 + 0(1). However, since 
T) G A(0), we have 77I = 0(p 2/3 ). Thus, 

* = (p 2 + O^ 3 )) 1 / 2 = p(l + 0(p" 4 / 3 ))^ = P (i + 0(p- 4 / 3 )) = p + O^ 1 / 3 )), 

which finishes the proof (recall that r\\ is positive). □ 

Since points in S 2 have the same first coordinate as the points from S 1; we immediately 
obtain: 
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Corollary 5.6. Let £ G H 2 (0). Then |^ - p| < p" 1 / 3 . 
Let us denote 

p_ = p _(0) : = inf{77i, ?7 = fa, r) 2 ) G Si(0)} 

and 

= P+(°) := sup-fa, 7? = fa, e Si(0)}. 
Then lemma 15.51 implies p+ — P- <C p -1 / 3 . Moreover, we can give another equivalent 
definition of S 2 : 

(5.12) H 2 (0) = {rj = (771,772), ?7i G (p_,p+)}. 
Note that we obviously have the following equalities: 

(5.13) p- = inffa, 77 = fa, 772) G H 4 (0)} 
and 

(5.14) a = inf{|p 2 |, 77 = fa,7? 2 ) G H 4 (0)}. 
Denote 

(5.15) p_ := supfa, 77 = fa, 772) G H 4 (0)} 
and 

(5.16) a := sup{|?7 2 |, 77 = fa, p 2 ) G H 4 (0)}. 
Lemma 5.7. W^e /iave: j5_ — p_ = 0(p~ l ) and a — a = 0(p^ 1 ^ 3 ). 

Proof. Let 77 := (p_, a), 77 := (p_, a), and 77 := a). Then all these points belong to 
H 4 (0). Thus, |f)| 2 — \ f]\ 2 = pt—P- = 0(1). Since p_ ~ p, this implies /3_ — p_ = 0(p _1 ). 
Similarly, |f)| 2 — |r)| 2 = a 2 — a 2 = 0(1). Since a ~ p 1 / 3 , this implies a — a = 0(p -1 / 3 ). □ 

Lemma 5.8. Suppose R n < pi/ 10 , £ G (5 3 (0) U 5 4 (0)), and let 7 G ©' 15A ~ f fre linearly 
independent of 0. Put 77 := £ + 7. T/ien ||t7| 2 — p 2 | ^> p 4//5 and, m particular, 77 ^ A(p). 

Proof. Since 7 and are linearly independent, 71 7^ 0; moreover, corollary 15.31 implies 
|7i| R~ 2 p -1 / 5 . Corollary 15.61 implies |£i — p| p -1 / 3 . Thus, |?7i — p| ^> p -1 / 5 and 

l^-p 2 |»p- 1/5 fa+P)»P 4/5 - 
Since ?7 2 <C p 2 ^ 3 , this implies ||?7| 2 — p 2 | 3> p 4 ^ 5 , so 77 ^ A(p). □ 

Now we make one more definition 

(5.17) H 5 (0) := H 3 (0) \ (H 4 (0) + ^ez{j6}). 

Lemma 5.9. Suppose £ G H 5 (0) and j G Z. If£ + j0 G A(0), i/ien £ + j0 G H 5 (0). 
Proof. Indeed, our assumptions imply that $, + jO G 2 3 (0). Moreover, 

^ +J 0^(H 4 (0) + U jez {j0}). 

□ 
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Lemmas 15.81 and 15.91 immediately imply 

Lemma 5.10. Suppose R n < pi/ 10 , £ G H 5 (0), and 7 G ®' 15J ^- If £ + 7 G ^L(p), i/ien 
| + 7 GS 5 (0). 

Proof. If 7 is linearly independent from 0, this is proved in lemma 15.81 Suppose that 
7 = j6, j G Z. Then £ + 7 G 2 2 (0), so if we assume £ + 7 G .A(p), this means that 
either £ + 7 G 2 3 (0), or £ + 7 G S 4 (0). The last possibility contradicts the definition of 
H5(0). Thus, £ + 7 G A(0) and now the statement follows from lemma 15791 □ 

Lemma 5.11. S x (0) C H 5 (0). 

Proof. Definitions of the sets Hj immediately imply that Si(0) C S 3 (0). Thus, it remains 
to prove that if £ G H 4 (0) and j £ Z, j ^ 0, we have 77 := £ + j0 G" A Without loss 
of generality we may assume that £2 > 0. Then £ 2 G [a, a] (see (15.141) and (15.161) ). and 
thus lemma [5T71 implies £2 — a + 0(p~ 1 ^ 3 ). The second condition in (I5.6P implies that the 
distance between each point of the set {a+j\9\,j G Z, j ' ^ 0} and ±a is at least Thus, 

||%| — IC2II > 3 f° r sufficiently large p. Since 771 = £1, this implies | |t7| 2 — |£| 2 | ^> p 1 / 3 ^. 
Since £ G A, this means that r) ^ A. This finishes the proof. □ 

Now we discuss the relationship between S 5 (0j) for various j. 

Lemma 5.12. Suppose, R n < pl/ w and ji ^ j 2 . Then (£5(0^ ) + 6 15J ^)nH 5 (0 j2 ) = 0. 

Proof. Denote rjj := pn(0j~). Then the definition of S5 and corollary 15.61 imply that the 
distance between any point £ G H 5 (0j) and 77^ is 0(p 1//3 ). On the other hand, lemma I5"72l 
implies that |j7 - a — ry - 2 | 3> -R^ 2 p 3> p 1 / 2 . Thus, if $ >1 G £5(0^) and £ 2 G 5 5 (0j 2 ), we have 
l£i — €2 1 ^ P 1 ^ 2 ■ Since p 1//2 ^> i? n , this finishes the proof. □ 

Now we define 

(5.18) V = <D(p) : = Uf =1 E 5 (0,) 
and 

(5.19) B = S(p):=yt(p)\0)(p). 

The sets H 5 (0) are called resonance regions corresponding to 0. The set D is called 
the resonance region. Finally, the set 23 is called the non-resonance region. Obviously, 
23 consists of L connected components, each one is located 'between' 0j~ and 0j~ +1 for 
some I, where of course we use the convention that 6l+i = 0i- We call this connected 
component (located 'between' Of and J+1 ) 23;. More precisely, we define (see figure 4 
at the beginning of the section 6) 

(5.20) 23, : = {x G 23, (£, 0,) > & (£, l+1 ) < 0}. 
We also define 

(5.21) ~ o (0) :=S 5 (0) + 6 7 ,- / 



DENSITY OF STATES 



23 



and 

(5.22) S (S) :=3 + 0^. 
Lemma 5.13. We have: 

(5.23) (20(^1) + e A )nSo(^) = 

when ji 7^ ]2 and 

(5.24) (So(S) + e^) n 5 (^) = 0. 

Proof. Formula (15.231) follows from Lemma [5. 121 Suppose formula ( 15.24ft does not hold. 
Then there exists a point £ G H 5 (0j) and 7 G @' 15A ^ such that r] := $, + j E H> C A. 
Lemma 15.101 implies that 77 G S 5 (0j). This means that 77 G" S in view of (I5.18P and 

(JEH- □ 

Let us introduce more notation. Let C C W 1 be a measurable set. We denote by 
J'W(C) the orthogonal projection in f) = L 2 ([0, 27r] d ) onto the subspace spanned by the 
exponentials e^(x), £ G C, {£} = k. 

Lemma 5.14. For arbitrary set C C M. d and arbitrary k we have: 

(5.25) v'? (k) (e) = y (k) (e + e 1 )^'? (k) (e) 

Proof. This follows from the obvious observation that if£ = m + kGC and |n| < R n , 
then £ + n G (C + Gi). □ 

We are going to apply lemma IBTTl and now we will specify what are the projections Pj. 
The construction will be the same for all values of quasi-momenta, so often we will skip k 
from the superscripts. We denote P l := 9^(E o (0i)), I = 1, . . . , L and P° := ?( k )(H (£)). 
We also put 

pj ■.= o>« ((Sb(o,) + e 6A?+ .) \ (s 5 (0/) + e 6 ^ +J ^)) , 1 = 1, . . . , L, j = 1, . . . , M, 

P^:=T( k )(S 5 (^) + e 6Ai ),/ = l,...,L, 
po ;= y(k) ^3 + e , } \ (s + e ; 7 = i ; . . . ; m ; 

P ° := 3>W(S). 

Finally, we define Q := I - Ef=o P '> #( k ) := Ef=o P ^'( k ) P ' + <2#o(k)Q, and J := 
[A-90v,A + 90v]. 

Lemma 5.15. Le£ /i„(P'(k)) G J. T/ien |/i„(P(k)) - fi n (H'(k))\ < p- 2 */ 3 = p- 2M . 

Proof. This follows from lemma 15.11 and from properties of the sets H formulated in 
lemmas [5 . 4H5 .111 Indeed, let us check that all the assumptions of lemma loTTl are satisfied. 
Lemma EH implies that if n then P l P n = P l V'P n = 0. The properties PjVPj = 
for \j — t\ > 1 and PjVQ = if j < follow from lemma 15.141 The distance from 
the spectrum of QH Q to J is greater than Av: this follows from the fact that Q is a 



24 



L. PARNOVSKI & R. SHTERENBERG 



projection to all the exponentials with £ lying outside of the union S (B) U UiE o (0i) 
and, thus, satisfying £ G" A. Finally, let us show that the distance from the spectrum 
of P!jH Pj (j > 1) to J is greater than cp}/ 3 . When / = 0, this follows from Lemma [5.41 
and the fact that A(0 P ) fl S = for any primitive vector p from ©l^- Suppose, I ^ 0. 
It is enough to prove that if 

77 G ((36(00 + 7 m) \ (Hb(0i) + %m)) 

with 6 1 G 6 gA ^, then 

(5.26) ||r/| 2 -p 2 | >p 1/3 . 

Since 77 G (5 5 (0/) + O 7J q-), we can write it as 77 = £+7 with £ G S 5 (0;) and 7 G 7A /. If 7 
and 0; are linearly independent, ( 15.261) follows from lemma [5781 Suppose, 7 is a multiple 
of 0\. Let us introduce coordinates generated by 0\ as above (after corollary [573]). Then we 
have jji = (1 6 [p_,p+]. Moreover, since 77 G" S 5 (0;)+©6M> we have 1 772 1 > a+|0j|. Indeed, 
suppose 1 772 1 < a + \0i\. Assume as we can without loss of generality that r] 2 > 0. Then 
77 - Oi G Sg(0,). Since £ G S 5 (0,), we have: 77 - ; = £ + G" 3 4 (0/) + Z0j. Therefore, 
r) — 0i E E 5 (0i), so 77 G 3 5 (0;) + 6A ^-. This contradiction shows that 1 772 1 > a + |0;|. 

Let us now denote by u the point with coordinates v x = p- and z/ 2 = a. Then 1/ G A, 
so ||i/| 2 - p 2 | < I. But 

It/I 2 - H 2 > ?? 2 - a 2 > (a + |0,|) 2 - a 2 > a > p 1/3 . 

Thus, 1 77 1 2 — p 2 ^> p 1 / 3 , which finishes the proof. □ 

Now we are going to construct mappings /, g : M 2 — > R with properties stated in the 
previous section. Let ^ G M 2 with {£} = k. Then we are going to define 

(5.27) /(£) = /i p (i/(k)) 
and 

(5.28) g(£)=p p (H(k)), 

where p = p(£) is a natural number chosen in a certain canonical way so that the mapping 
p : G 1R 2 , {£} = k} — > N is a bijection. Leaving aside for a moment the question of 
the precise definition of this mapping, we notice that if we define the functions / and g 
by formulas (15.271) and (I5.28p . then the properties (i) and (ii) formulated in the previous 
section will be satisfied due to lemmas 15.151 and 15.11 So, now we discuss how to define 
the mapping p. Before doing it, we need more definitions. Let £ G A. Then £ belongs 
to exactly one of the sets S, 3 5 (0i),..., E 5 (0 L ). If £ G 25, we define 

(5.29) T(€):=* + e A . 
If £ G E 5 (0i), we define 



(5.30) 



r(0 = t(£ 00 := + jo, g 3 3 (0,), iez} + e 7A ^ 
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We call two vectors £ x and £ 2 equivalent, if T(^ x ) = Y(£ 2 )- Note that £ 1? £ 2 e ^ could 
be equivalent only if they belong to the same H 5 (0j). 

Now suppose that 77 G M 2 . Then we can define T(t7) in the following way: if 77 G Ho (23), 
then we have T7 G Y(£) for a unique £ G .A (then £ G 23); if t] G H (#z), then we have 
77 G Y(£) for a unique (up to the equivalence) £ G D (then £ G H 5 (0/)). In both these 
cases we put Y(r?) := Y(£). Finally if 77 G" (S (S) U Uf =1 H o (0;))> we put Y(r/) := {rj}. 
We also define P(rj) := 0>«»»»(T(f7)). 

Lemma E3] implies that the operator P°P'(k)P° admits a decomposition into invariant 
subspaces: 

(5.31) P°H'(k)P° = P(£)H'(k)P(£). 

£e£,{£}=k 

Similarly, lemma [5T81 implies that for each / = 1, . . . , L we have: 

(5.32) P l H'(k)P l = P(£)Pf'(k)P(£), 

where the union in the RHS is over all classes of equivalence of £ G Hs(0j) with {£} = k. 
Finally, we obviously have: 

(5.33) QH (k)Q = P(rj)H (k)P(ri), 

where the union is over all 77 G" (H (23) U Uf =1 H Q (0i)), {77} = k. Moreover, since all 
projections P(r7) in (15.331) are one- dimensional and we have assumed that f Q V(x)dx. = 0, 
we can replace P^o(k) with H'(k) for the sake of uniformity so that 

(5.34) QP (k)Q = P{r,)H'Qt)P(<n). 
Thus, 

(5.35) P(k) = 0P(r7)P'(k)P(r7), 

where the union is over all (non-equivalent) 77 G M 2 , {77} = k. 

Suppose now 77 G IR 2 , {77} = k. Then |?7| 2 is an eigenvalue of P(/7)Po(k)P(r7), say 

(5.36) |r7| 2 = ^(P(r 7 )P ! o(k)P(r7)). 

If \r]\ 2 is a simple eigenvalue of P(rj)H (k.)P(rj), then this defines the number t uniquely. 
Suppose now that |r7| 2 is a multiple eigenvalue, say \t]\ 2 = \i)\ 2 , 77 G Y(ry), and there 
are precisely t — 1 eigenvalues of P(r])Ho(k)P('q) below |t7| 2 . In this case, we label these 
eigenvalues according to the crystallographic order of their universal coordinates. More 
precisely, we write \r)\ 2 = fj lt (P(7])H (k)P(rf)) and |t)| 2 = fi t+ i(P(rj)H (k)P(r])) if either 
r)i < fjx, or rji = fji and r] 2 < fj 2 . Thus, we have put into correspondence to any point 17 
a number t = t(rj), t varies between 1 and the number of elements in T(rf). (Although 
we will not use this function t(rj) in this section, it will be of much use for us later on). 
Next, we define 

(5.37) v^) := ^{P^H'^P^)). 
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Due to (15.35p . the set {v(t]), {t}} = k} coincides with the set of all eigenvalues of 
H(k) (including multiplicities). Let us label these eigenvalues in an increasing order; 
in the case of multiple eigenvalues we, as before, label them in accordance with the 
crystallographic order of their coordinates. Then to each point r], {77} = k, we have put 
into correspondence a number p = p(rj) such that 

(5.38) u{ri) =fi p (H(k)). 

Thus defined mapping p is the mapping we are using in the definitions (I5.27P and (I5.28P . 
The rest of this paper is devoted to introducing the coordinates (r, $) and checking that 
the conditions of lemma H. 101 are satisfied. We start from the non- resonance region S. 

6. Non-resonance regions 
Suppose that £ G 23/ (recall that 23/ is defined in ( 15.201) and </>(x 1 ,x 2 ) is the angle 



between two non- zero vectors and x 2 ). Put 



(0i,fli+i) 



Throughout this section, 



we fix the coordinate (£1, £2) introduced after corollary 15.31 and related to Q\. Namely, we 
put £1 = (£, n(0/~)) and £2 = (£) n (0/))- There is a unique point v = v{l) satisfying the 
following two properties: (u,n(0i)) = a and (u, n(0i + i)) = —a; we have z/i = acot0z, 



l>2 



a, so \v\ 



sin 4>i ' 




Figure 4 

We introduce the following pseudo-polar coordinates (r, <E>) on 2^: r(£) := |£ — ui\ 
and $(£) = 0(£ - u,0j-) when £ G Obviously, G [0, 2&] =: §/ = Sp and 

r(£) ~ p when £ G 23;. We also have the following formulas: £1 = + r(£) cos($(£)) 
and £2 = ^2 + KO s i n ( < ^ ) (^))- Therefore, 



(6.1) 



(r(£) cos($(£)) + \v\ cos fa) 2 + (r(£) sin($(£)) + Hsin<^) 2 

(r(£) cos($(£)) + a^) 2 + (r(£) sin(<&(£)) + a) 2 , 
sm </>, 
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This implies that there is a complete asymptotic formula: 

/ oo 

(6.2) 1 + 



3=1 



~ ■ i 

with bj = bj <C a 3 Rn ^ pn as r(£) — > oo, uniformly over £ G !Bj, and this formula 

1 /12 

can be differentiated once. (Here we assumed that R n <C p n •) 
The following lemma was proved in [8] (Lemma 6.1 there): 

Lemma 6.1. Let R n <C pi/ 24 . Then the following asymptotic formula holds: 

<?(£) ~ \t\ 2 



(6.3) 



+ E E E m.&ih> mi ■■■&»!.> 

s=l r]i, -,V s £®'-, m 1 + -+m s >2 



in a sense that for each natural K we have 

+ E E E A m 1 ,...,m s (cvi)- mi ---(tv s r ms +o(p- K ) 

s=l t?, ,...,T7„e0'- miH \-m B >2 

uniformly over R n < Pn /24 and £ el #ere ; A miv . . m is a polynomial of the Fourier 
coefficients V(rjj) and V(f]j — f]i) of the potential and the exponents mi, . . . ,m s are 
positive integers. Moreover, 

(6-5) \Am,...,m p \ < 1 

uniformly over n (but with the implied constant depending on V and m\,...,m v ). 

Remark 6.2. Estimate (16.51) was not stated in [8], but it follows easily from the proof 
of Lemma 6.1 there. 

Corollary 6.3. We have: 

(6.6) ^) = l€r + G(€) + o(p- Jf ), 
where 

M 

(6.7) g«):=e E E ^..^(cmr™ 1 ---it, ri s r m °. 

s=l ■n 1 ,...,7} s eQ' & 2<m 1 +---+m s <M 

Remark 6.4. Since, as we have seen in lemma H~3l and corollary 14.41 the terms of order 
0(p~ M ) do not contribute to asymptotic formula (13.61) . we can re-define 

(6.8) g(t):=\t\ 2 + G{Q. 
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Lemma 6.5. Assume R n <C pi/ 20 . For each in G N and r\ G such that rj is not a 
multiple of Oi or Oi + i there is a complete asymptotic formula: 

oo 

(6.9) (£r7>-™~£r(£)-^(<&(£)) 

j=m 

uniformly over £ G "Bi, where \df l \ -C pi! 2 . Similar formulas are valid if rj is a multiple 
°f an d < $(£) < 0/, or if r] is a multiple of Oi and 0; < $(£) < 20/. 

Proof. We have: £ = f + r(£)n(£ — i/). Therefore, (£, 77) = (1/, 77) + r(£)(n(£ — v),rj). 
Our constructions and corollary 15.31 imply |(n(£ — t-') , r/) j ^> i?~ 2 . Thus, 

(£, 77)- 1 = r-^riH - u), n)-\l + r-^riH - u), r,)~ l {u, r,))- 1 

(6 ' 10) =r- 1 (n(£- l ,),r 7 )- 1 f;(-l)^(n(£- l ,),r 7 )^(^,r7V 

i=o 

and (n(£ — u),rj)~^(u,rjy <C pi/ 3 <C p^ 2 . Now formula (I6.9P is obtained from 
(16.101) by raising both sides to the m-th power. The proof of the last two statements is 
similar. □ 

Unfortunately, lemma 16.51 does no longer hold if 77 is a multiple of Q\ or 6i+\ and 
$(£) is close to or 20/ respectively. Therefore, we cannot apply lemma 14.101 without 
modifications. This means, we need to do some extra work. We can assume, without 
loss of generality, that 0; < 1/100, which is certainly the case for sufficiently large n. 

Let us fix an angle $ for a moment, < $ < 20; and let £ = (r, $) where only r 
varies. Denote by r = r ($;p) a unique value of r which corresponds to £ satisfying 
|£| 2 = p 2 . It is easy to check that the partial derivative |^ = 0(p~ 4 / 3 ). Therefore, there 
is a unique value of r such that corresponding point £ = (r, $) satisfies </(£) = p 2 ; we 
denote this value of r by r\ = ri($; p). 

Lemma 6.6. There is an asymptotic decomposition 

00 

(6.11) r ~p(l + ^p^), 

3=1 

where pj = Pj(&) = 0(pi! 2 ) uniformly over $. Moreover, we have p\ = — acos(0/ — 
$)(sin0;)- 1 ; p 2m = (^ 2 )(-l) m a 2m sin 2m (0, - $)(sin 0/)" 2 "\ and p 2m+1 = for m G N. 

Proof. This follows from the explicit formula which can be easily obtained using the 
cosine theorem: 



(6.12) r = -acos(0; - $)(sin0;) 1 + J p 2 - a 2 sin 2 (0/ - $)(sin ^~ 2 



•i) 

□ 
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Lemma 6.7. We have: 

i r*k 

(6.13) vol(i+nS z )-vol(i-n30 = - / (r 1 ($;p) 2 -r ($;p) 2 ) C /$. 

^ Jo 

Proof. Integrating in polar coordinates, we have: 
(6.14) 

/■2</>i rri r2<j>i pro r2<pi 2 _ 2 

vol(i + n B z ) - vol(i _ n S { ) = / d$ / rrfr - / d$ / rrfr = / ^ 

Jo Jo Jo Jo Jo 2 

□ 

The last two lemmas show that in order to compute vol(A + PI Hi) — vol(v4~ H 2;), it 
remains to compute r\. We do it using the sequence of approximations. Assume as above 
that $ is fixed. Put tq := r$ and £ := (fo, The further elements of the sequence are 
defined like this: £ m+1 = (f m +i, is a unique point satisfying |£ m+1 | 2 = P 2 — G($, rn ). 

Lemma 6.8. For each m e N we /iaue: 

(6.15) r 1= f m + 0(p- m ). 

Proof. Put 



(6.16) #(r) := -acos(0; - $)(sin<^) _1 + \J p 2 - G(r, $) - a 2 sin 2 (0, - $)(sin^)- 2 . 

Then H'(r) = 0(p~ 4 / 3 ). Moreover, n is a unique solution of equation ri = H(r\). 
Thus, Banach contraction mapping theorem tells us that the sequence r m satisfying 
f m+ i = H(r m ) converges to n and |r*i — f m+ i| <C p _1 |?"i —r m \. Since ri = r + O(l), this 
finishes the proof. □ 

Corollary 6.9. We have: 

(6.17) vol(A+nS,)-vol(A-nS,) = r / (r M+1 ($;p) 2 -r ($;p) 2 )d<l> + O(p- M ). 

^ Jo 

Analogously to (16.12j) . we have: 



(6.18) r m+1 = -acos(<fo - $)(sin<^) 1 + yj p 2 - G(r m , $) - a 2 sin 2 (0, - $)(sin0,) 
Taking into account (16. 7p . (16.101) . (I6.12p . (16.181) . and lemma [631 we obtain that 



r M +i = r 



(6.19) " J " 

= r + P" 1 £ C jAM (Q)p->(S , n(^))" s + 0(p 



j,s>0;2<j+s<2M 
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for < $ < 4>i. Similarly, 

fM+^ro + p- 1 Yl A>,M($)p^(| >n(^+i))" s 

j,s>0;j+s>2 
6.20 " " 

= r + p- 1 Yl C^M^P-Hto, n(0,)>- + O^- 1 ) 

j,s>0;2<j+s<2M 

for (pi < $ < 20/. Here, Cj iS) m($) and C^mO^) are polynomials of cos(<p; — $), sin(<p; — 
$), and expressions of the form ((n(£ — is), n(T/ t ))) _1 , where r] t G are not multiples 

of 6 1 or 9i + i. Each term in the polynomial Cj )S ^i{^) and C/ )a ,Af($) is 0(pn 2 ) } and the 
number of such terms is 0(Rn^ +s ^)- 

Next, we note that (£ o ,n(0/)) =a + r sin$and (£ o ,n(0/ +1 )) = — a — r sin(20; — $). 
Thus, in order to use Corollary 16.91 we need to compute integrals of the form 

(621) [ C^m^M*) 6 d* f C jia , M (*)ro(*) 6 d* 

1 ' J 7 (a + r ($)sin($))*' 7 (-a - r ($) sin(2^ - ' 

o w 

Taking into account Lemma 16.61 and properties of (7j iS) m(3>), C^>,m($) stated above, 
we can decompose all functions of $ in (I6.2ip into Taylor's series in the neighborhoods 
of $ = and $ = 2<p/). At the same time we apply the following transform of the 
denominator: 

(a + r ($) sin($))- s = (a + p$(l + p(p, $))(1 + </)($)))^ s = 
/ $ 

(6 22) (a + P$) " V + a/p+$ (p(P ' $) + + P(P ' 



where (see Lemma 16. 6ft 

oo 

p(p,®) ■= Y pk ^p' k 

k=l 

and 

= Y" 1 l) $ 2fc - 
^ (2^ + 1)!^ ' 

recall that we are assuming that (pi < 1/100, so that there is no doubt about the conver- 
gence of the last series in ( 16. 22ft . Thus, decomposing pk{<&) into Taylor's series we reduce 
the problem to computing the following model integrals: 

{a/p+ <£>) m ' 
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After substitution x := a/p + $ we can explicitly calculate these integrals. Note, that 
ifl<m<A; + l then the term In p appears. Combining together all contributions we 
obtain the following lemma, which is the main result of this section: 

Lemma 6.10. Assume that R n <C pll 2A . Then 

M M 

(6.23) voi(i + n - voi(i- n %) = J2 c iP~ 3 + ln P Qp~ j + °(p~ m ), 

3=1 3=2 

where \C X \ < pn l/ \ |C,-|, \CA < p 2 n j/3 , j > 2. 

7. Resonance regions 

We now consider 77 G D and try to compute g{rf). The key result in this section is 
corollary 17.111 where we compute g{rj) in this setting. In the rest of this section, we 
fix I and omit it from the notation, so that := B\. We also assume that n is fixed 
and will frequently omit it from the notation. As above, we introduce the coordinates 
77 = (771, 772) so that rji = (r],ii(6 ± )) and 772 = (77,11(0)). We obviously have 771 ~ p and 
1 772! p 1 / 3 . It is convenient to denote r := 771 and $ := 772, to indicate that (r, $) are 
going to play the same role as in corollary 14.41 (or rather remark 14.61) and lemma 14. lOt 
note that (r, <3>) satisfy all properties of remark H~6l We also fix an element £ G Hs(0/) 
in each set T and assume that 77 G Y(£); the point is that we will frequently treat £ as 
fixed and study how g{rj) varies when 77 runs over Y(£). 

Let v>q = 0,i/i, ... ,u p be a complete system of representatives of 9 7J ^ modulo 0. 
That means that Uj G 7A ; f and each vector 7 G 9 7J q- has a unique representation 
7 = Vj + m0, m G Z. We denote the coordinates of z/j by (z/,z/') and put Vtr, = 
:= (£ + 17 j + (^)) l"l Y(£). Then each set consists of points having the same 
first coordinate; the distances between points in SE^ are multiples of |0|. Moreover, 

(7.1) T(€) = U*i, 

3 

and this is a disjoint union. 

Let us compute diagonal elements of H(£) := P(£)H'(k)P(£), where k = {£}. Put 

:= P(£)f), so that H(£) can be thought of as an operator acting in $)(£). 
Let 77 G Y(£). Then 77 can be uniquely decomposed as 

(7.2) 77 = £ + m0 + uj 

with m G Z. Recall that H(£) = P(£)(H Q (k) +V')P(£) and H ($,)e v = \rj\ 2 e v whenever 
77 G Y(£). We obviously have: 

|7y| 2 = |£ + Uj + m0\ 2 = {r + v\) 2 + + i/j + m\0\) 2 

= r 2 + 2u' f + vf + (6 + < + m|0|) 2 - 
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This simple computation implies that 

(7.4) H(£) =r 2 I + rA + B. 

Here, A = A(£) and B = B(£) are self-adjoint operators acting in P(£)j} in the following 

way: 

(7.5) A = 2j2" , j y ( - k) (* j ); 

3=0 

in other words, for r\ G *S?j we have 

(7.6) Ae v = 2^e„ = 2(rj - $)ie v , 
and 

(7.7) Be v = {vf + (6 + ^ + m\6\f + P(0^'K 

for all r\ G ^j(^) with v 3 and m being defined by (17.21) . These definitions imply that 

(7.8) QJ :=ker A = ?< k > (¥(,)£(£)• 
Notice that 

(7.9) P(|)||« J R n « P y 3 
and 

(7.10) ||P(£)|| «P^ /3 . 
Let us state more properties of A and P. 

Lemma 7.1. Le£ /j, be a non-zero eigenvalue of A. Then 3> P^ 2 - 

Proof. Formula (17. 6p implies that the eigenvalues of A equal {Uj}- We also have: z/j = 
(uj,n(0 )). Now the statement follows from corollary 15.31 □ 

Let us define P to be the orthogonal projection onto 23 = kerA acting in S)(£) and 
B := PBP : 23 — > 23. Note that considering operators acting in 23 means considering 
only j = (and thus Uq = 0) in ( 17.51) and ( IT. TP . Thus, in particular, we have: 

(7.11) Be v = ((6 + m\0\f + P(£) V')e„ 

if 77 = $, + mO G \&o (^) • We also denote fi2 := and hi '■= {rar } (note that n2 is not 
the second coordinate of n = [£]; this is why we did not call it n 2 ). 

Lemma 7.2. We have: 

(7.12) // i+2 (P)-^(P)>l 
uniformly over j , n, I and £ G H 5 (0). 
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Proof. Denote by T the number of elements in {% + jO, j G Z} H A(0). Inequality ( 17.12ft 
obviously holds if j > T — 2. Indeed, denote by L? the operator B with potential V 
being identical zero. Then we have \fij[B) — fj,j(Bo) \ < v. On the other hand, it is easy 
to check that yUj +2 (_B ) — ^j(Bo) ^ a ~ P 1 ^ 3 - 

Let us assume now that j < T — 2. Then we will compare eigenvalues of operator B 
with the eigenvalues of a certain one-dimensional Sturm-Liouville operator. Let Y(£) := 
{€+j0, j G Z}, P(£) := ?( k )(f (£)), and := P(£)#- Consider an operator B = B n 
(later on in the proof, we will need to remember that these operators depend on n) acting 
in by the formula 

(7.13) Be v = ((( 2 + j\0\) 2 + P(t)V')e r , 

for each 77 = £ + jO G Y(£). Then, in the same way as we proved lemma 15.151 using 
lemma EUl we can show that if j < T, we have 

(7.14) \ N {B)- N {B)\ «^- 1 V 3 . 

However, the operator B is unitary equivalent to a one-dimensional Schrodinger operator 
— y" + V on the interval [0, 27r|0|~ 1 ] with a potential 

(7.15) V = V e , Rn = Yl 0^) V2 e lx ^ +mW) V(m6) 

m£Z, \m0\<R n ^ ' 

and quasi-periodic boundary conditions y(j^) = e 27Tk2l y(0) and y'(^r) = e 2nk2l y'(0). 
Indeed, the isometry S which establishes this unitary equivalence is given by S : e^ +m e 1— > 

1 /2 

m. j gia(6+w.| l). Standard results about one-dimensional Schrodinger operators (see 



2tt 

e.g. 1X0] ) imply that 

(7.16) N+2 (B)- N (B)->\6\ 2 . 

The simplest way to see why this inequality holds is to notice that the distance between 
eigenvalues of B and the unperturbed eigenvalues {(m + k,2) 2 \0\ 2 } m& i is at most the 
Loo-norm of the potential Ve,n n - This shows that (17.161) holds when j > Cv, whereas 
for finitely many j satisfying j < Cv we can use the fact that fij + 2(B) ^ H>j(B), since 
an eigenvalue of a one-dimensional differential operator of second order cannot have 
multiplicity three. Inequalities (I7.16P and (17.141) prove (17.121) for j < T. Let us prove 
that this estimate is uniform in j, n, and I. Indeed, the uniformity of (17.141) follows from 
lemma 15.11 Consider ( 17.161) . Uniformity in j follows from the remark after ( 17.161) . It 
follows immediately from the same remark that (17.161) is uniform when Loo-norm of the 
potential Ve,R n satisfies 

(7-17) ll^idloo 
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Since the potential V is infinitely smooth, we have | V r ( , y) | <C |t|~ 2 , which shows that 
there are only finitely many 6 for which ( 17.17P is not satisfied. This shows uniformity of 
(17.161) in I. It remains to prove the uniformity of ( 17. 16ft in n when is fixed. First, we 
notice that (17.161) holds for sufficiently large j > j , where jo depends only on 1 1 V| |oo, 
but not on n. Suppose now that (I7.16P is not uniform in n. Then there is a value of j 
such that 

(7.18) lim n j+2 (B n ) = lim fi j+1 (B n ) = lim fij(B n ) =: \i 

n— too n— >oo ri— >oo 

(strictly speaking, we need to pass to a subsequence n k if necessary). However, these 
limits are the eigenvalues of the limit operator with the potential 

(7.19) V e ,oo = J2 (|J) ' e lx ^ +mW) V(mO). 

The required result now follows from the fact we already used above that a second 
order one- dimensional differential operator B^ cannot have an eigenvalue of multiplicity 
three. □ 

Our next task is to compare eigenvalues of H(£) and H(£) when £ and £ are two 
different vectors lying in S 5 (0). This is not a straightforward task, since these operators 
act in different Hilbert spaces (S}(£) and £)(£') correspondingly). Thus, first of all we 
need to be able to map these Hilbert spaces onto each other. The natural candidate for 
such a mapping is 

(7.20) Ftf{e n )=e n rt_e. 

Ideally, we would like this mapping to act as follows: F^g : fj(£) — > -£)(£') an d be an 
isomorphism. Unfortunately, in general this is not the case since the sets Y(£) and Y(£') 
can contain different number of elements. In fact, it may well happen that rj & Y(£), but 
(77 + £' — £) ^ Y(£'). However, the mapping F has the suggested property in one very 
important special case: when = (in other words, when the second coordinates 

of £ and £' coincide). Indeed, suppose that $(£) = Then obviously 

{t + jOt e S 3 (0,), j'eZ}+ (£' - £) = tf + jOt e 3 3 (0,), j e Z}. 
Thus, we also have 

T(0 + (e'-0 = T(O, 
and so the mapping is an isometry between $)(£) and ■$}(£') with = Fg£. 
Moreover, if we look carefully on formulas (I7.6P (the first equality there) and (17.71) . we 
realize that the definitions of operators A(£) and -B(£) do not depend on £1, so we 
have Fg £A(£)Fgg = A(£) and, similarly, Fg ^B(^')F^ g = B(£). Thus, all operators 
A(£) are unitary equivalent when £ runs along any horizontal line = the 

same statement holds for B(£). It is convenient to think of all such operators as being 
identical operators A($q) and -B($ ) acting in the same Hilbert space fj($o). We also 
notice that if $(£) = $(£')> then the isometry Fg ^ leaves the function t (defined after 
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( 15.361) ) invariant. This means that whenever 77 e Y(£) and 77' = 77 + £' — £ e Y(£ ), 
we have £(77) = £(77'); this is true even if |r/| 2 is a multiple eigenvalue of H (r)) (and, 
correspondingly, |t7'| 2 is a multiple eigenvalue of H (rj')). 

Denote § = S^ +L := [—a, a], / = 1,...,L (recall that §' were already introduced in 
the previous section). Now it seems to be a straightforward task to apply lemma 14.101 
in the resonance region similarly to how we did it in the non-resonance region. Indeed, 
formulas (I4.9j) . (I4.10p . and (14.121) are immediate corollaries of |£| 2 = r(£) 2 + $(£) 2 , and 
(14.111) follows from the standard results of perturbation theory (see, e.g., [5]) applied to 
the operator pencil rA(<$>)+B(<&) = r(A($)+r _1 i?($)). The problem with this approach 
is that the coefficients Oj($) in (14. lip are not bounded in general. This unboundedness 
of the coefficients is caused by the fact that the eigenvalues of B can be located very 
close to each other. However, lemma 17.21 shows that the multiplicity of any cluster of 
eigenvalues of B cannot be greater than 2. This observation will be of a great help to 
us. 

It will be slightly more convenient to introduce new operators A = A n := P? n A and 
B = B n := P? n B (and B := P? n B = PBP); we will be assuming from now on that 
R n < pi/ 25 . The reason for this change is that lemmas 17.11 and [ 7.21 can be reformulated 
in a more uniform way: 

Lemma 7.3. There is a positive constant C 2 which satisfies two properties: if fj, is a 
non-zero eigenvalue of A, then > C 2 and 

(7.21) fij +2 (B) — jUj(B) > 3C 2 

uniformly over j , n, I and 

Remark 7.4. 1) It will be convenient to assume that C 2 < 1/10, which we will be doing 
from now on. 

2) Of course, we have slightly better estimate for eigenvalues of /ij(B). The distance 
between yu J+2 (B) and /U,(B) is 3> P 2 . But (17.211) is enough for our purposes. 

The importance of lemma [7731 can be seen from the following remark. Suppose that we 
could establish the inequality (17.2ip with /i J+1 (B) instead of /ij +2 (B). Then, using the 
approach from the previous section, we could prove that the coefficients %(<!>) in (14. lip 
are bounded, and this would finish the proof of our main theorem. However, in general 
it could happen that two eigenvalues of B lie close to each other. Our further course of 
action will reflect this possibility. We will divide the segment S = [—a, a] into two disjoint 
parts, S = S U §. Roughly speaking, S will be the region where the eigenvalues of B 
are far from each other, and § will be the region corresponding to couples of eigenvalues 
of B lying close to each other. To be more precise, we need yet more notation. Let 
£ G S 5 (0) and 77= (771,772) G *o(£) (recall that * = {£ + Z0}nY(£), so P = ? (k) (*o); 
this means, in particular, that 771 = £1). Then (772) 2 is an eigenvalue of the unperturbed 
operator B (£), say (ffc) 2 = A i r(T 7 )(Bo). Here, as above, we use the convention that if two 
eigenvalues (rj 2 ) 2 and say (y 2 ) 2 coincide, we label them according to the crystallographic 
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order of their universal coordinates (of course, this could happen only if the quasi- 
momentum A; 2 is either or 1/2). Thus, we have defined a mapping r : \l/o(£) ~~ > N. 
Notice that we can talk simply about the value r(ry), without specifying what £ is, since 
if 77 G *o(£j), j = 1, 2, then *o(£i) = *o(£ 2 )- Next > for an y P oint *7 G *o(£)> we define 

(7.22) h(ri) := lirwiBit)). 
and 

(7.23) h(r 7 ):=^ T{T , ) (B(0) = i?^(r/). 

Then, we can reformulate ( 17.2ip like this: for each 77 G \^o(£), there is at most one point 
v G *()(£)) v 7^ V sucn that |h(i/) — h(r/)| < 3C 2 . Notice that this whole construction 
does not depend on the first coordinate T]i (we can recall the paragraph after the proof 
of lemma E2 at this stage), so we can think of r as a mapping r : 772 1— > t(t7), where 77 
is any point with second coordinate 772 such that 77 G ^o(£) f° r some £ G H 5 (0). Then 
the domain of thus defined mapping r is some interval / which consist of all second 
coordinates 772 of points 77 G ^o(^) with £ G H 5 (0); obviously, / D [—a, a]. Similarly, 
/i : 772 I— > /j r ( JJ8 )(B(^)) is a well-defined function on /. Let s be a small parameter which 
we will fix later on. At the moment, we put s = |C*2, but we will decrease s later. We 
define I\ = Ii(s) to consist of all points 772 from I such that there exists a non-zero 
integer m such that 7/2 + m\6\ G / and 

(7.24) j h(77 2 ) -h{r} 2 + m\0\) \< s. 

In other words, I\ consists of all points 772 such that the eigenvalue of B corresponding 
to r/ 2 is close to being multiple. We also put I = I (s) := I \ h(s). Let us study the 
properties of this partition. First of all, due to lemma I7T31 for each 772 G Ji, equation 
(17.241) is satisfied for precisely one value of m; obviously, then 7/2 + m\0\ also belongs to 
I\. Thus, we can define a mapping l : Ix — > Ji by the formula i(7/ 2 ) = ^2 + fn\0\, where 
m 7^ is chosen so that (I7.24p is satisfied. Obviously, then t 2 = Id. We can extend the 
mapping 1 to the whole / by requesting that ^(772) = 7/ 2 whenever 7? 2 G Iq. Sometimes we 
will slightly abuse this notation by writing i(r)) := (771, ^(^2))- 

Lemma 7.5. Suppose r/ 2 G [—a, a]. T/ien i(r/ 2 ) G [—a, a]. 

Proof. If r/2 G io, the statement is obvious. Suppose, 7/ 2 G (ii PI [—a, a]). Without loss 
of generality we can assume that 7/2 is positive. Notice that \h(r]2) — \i]2\ 2 \ < v. This 
implies that whenever 7/ 2 < a/2, the statement holds. Suppose, 7/ 2 > a/2 ^> p^ 3 . Then 
if (I7.24p is satisfied, we have 

(7.25) | \r] 2 \ 2 - 1 772 + m|0|| 2 |<s + 2?j, 
and thus 



(7.26) 
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Since 772 is assumed to be positive and rj 2 + m\0\ is negative (otherwise there is no chance 
for (I7.26P to hold), this means 

(7-27) I V 2 + ^\0\ |«p; 1/3 - 

Obviously, we will have the same inequality for ^.(772): 
(7-28) I \ L ( V2 )\ + ^\0\ |«p-V3 ) 

so both 772 and ^(772) are close to (i.e. within distance o(l)) points of the form ±y|0|. 
Now the second condition f!5.6p implies that the distance from a to any point of the form 
±yI#| is at least |0|/4. Thus, if \rj2\ < a, this implies that 1 1 (772 ) | < o,. D 

Now let us establish the relationship between the labeling r : ^o(^) ~~ * N used to 
define mapping h and the labeling t : Y(£) — > N defined by (15.361) . 

Lemma 7.6. Let T be the number of elements in Y(£) whose first coordinate is strictly 
less than £1. Then for each 77 G V&o(£) we have 

(7.29) t{rj) = T{rj)+f. 

Proof. It follows from the proof of lemma 15.81 that whenever \i G Y(£) \ ^o(C) an d 
77 G the following two conditions are equivalent: \i\ < r/i and < |?7|. Indeed, 

to say that fi G" \^o(£) is equivalent to saying that H\^r\\. Suppose that \i\ <r\\. then 
T]i — iii > i?~ 2 , so (?7i) 2 - (p x ) 2 > Pn-R^ 2 > Pn /5 . Since |p 2 | < a < pi*/ 3 , this implies 
|/Lt| < |r/|. The case pi > 771 is treated similarly. The rest follows from the definitions of 
mappings t and r. □ 

Now let us recall that because of (17. 4p . we are interested in studying eigenvalues of 
the operator pencil 

(7.30) Z(r) := rA + B = r(A + r _1 B) = i£#(r), 
where 

(7.31) Z(r)=rA + B, 

r = £ 1; and operators A and B depend on $(£) = £2 for some point £ = (£ 1; £ 2 ) £ 5 5 (0). 
These operators act in the Hilbert space which we have denoted by ($(£)); see the 
paragraph after the proof of lemma 17.21 for the discussion of this Hilbert space. To 
be more precise, we fix a point 77 G Y(£) and study the function girf) . We are only 
interested in eigenvalues of Z(r) which are bounded as r — > 00 (which means, they can 
be considered as perturbations of zero eigenvalues of A or, equivalently, that 77 G \l/o(£))- 
The first order of approximation of such eigenvalues as r — ► 00 are eigenvalues of B. We 
will consider separately two cases: 772 G Ji and 772 G Iq. The former case is much more 
difficult, and we give all necessary details. The latter case is much simpler and can be 
treated analogously to the first case (alternatively, one can apply methods similar to 
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those we used in the previous section); we will make some remarks on this case later. 
So, let us assume that r] 2 G I\. 

Recall that P is a projection onto ker A = ker A; we also denote P' := P(£) — P. Let 
Pq (< P) be projector onto span of two eigenfunctions of B corresponding to 11(772) and 
hW^)); put Pq := P—Pq. By fi we denote a spectral parameter, which at the moment we 
assume satisfies \fx\ < 2||B||. Operator P'AP' + P is invertible and IKP'AP' + P) -1 )! < c 
with constant c > uniform with respect to n and 772. Thus, for sufficiently large po, 
operator 

(7.32) 5 M := P'AP' + P + ^(B-fi) 
is invertible. We have 

(7.33) rA + B - fx = (I - PS' 1 )^. 

Then /i is an eigenvalue of rA + B if and only if is an eigenvalue of I — PS^ 1 . Obviously, 
corresponding eigenfunction 2 belongs to P5)(<3>). Thus, 

(7.34) 2 = S~ l z + y 
for some y G P'fj($). We have S^z - z = S^y. Then 

(7.35) P'S.Pz = P'S^P'y, y = {P' S .P'y 1 P' S ,P z . 
Thus, 

= (Ps,P -p- Ps.p'iP's.p'y'p's^z = 

(7.36) 1 _ _ 1 _ 1 

-(P(B - fi)P - -P(B - fx)P'{P'AP' + -P'(B - ^PO-'P^B - fi)P)z. 

Therefore, zero is an eigenvalue of operator P(B — fx)P — where 

(7.37) K„ = P(B - fi)P'(P'AP' + -P'(B - //)P') _1 P'(B - fi)P. 

r 

Note that since PBP' = R 2 n PV P' we have 

(7-38) \\K^\+j^\\f-K,\\<C{N)Rl 

where C(N) depends only on N and V, provided r > p^ 4 and po(V) is sufficiently large. 

Next, we want to narrow the range of /i's which serve as the candidates for being the 
eigenvalues of r A + B. Let us at the moment only look for eigenvalues fx such that 

(7.39) \fi-h(rf 2 )\<C 2 . 
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If we assume that ( 17. 39ft is satisfied, the operator P(B — p)P + P is invertible (on 03) 
and inverse operator is bounded uniformly in n and 772- We have 

1 k 

p, — v- 1 ~~ ± u-l^^ 



(7.40) P(B - ij)P - -K^ = (I- P,D- l )D 



where := P(B - p)P + P - ^ 

Now we repeat the same construction as in (17.34I) -( J7.37I) . only with respect to the pair 
of projections Po, P instead of P, P(£). As a result, we obtain that /1 is an eigenvalue of 
Z(r) = rA + B if and only if zero is an eigenvalue of the operator Po(B — p)Po — G M , 
where 

(7.41) G, = ^(PoAVPo + ^PqKMP&B - - ^P^P^P^K.Po). 

The operator Po(B — p)Po — G M is, in fact, a (2 x 2)-matrix. Note that in a suitable 
basis, Po(B — /i)Po is a diagonal matrix with h(r/ 2 ) — p and h(i(r/ 2 )) — P standing on the 
diagonal. Calculating the determinant of Po(B — p)Po — G M in this basis, we obtain that 
p satisfying (17.391) is an eigenvalue of Z(r) if and only if 

(7.42) (h(t&) - p + iai)(h(i(^)) - /x + -a 2 ) - ^/? 2 = 0- 

Here, ai, a 2 , and (3 are functions of p and r (depending on r/ 2 as a parameter) analytic 
in |p — 11(772)1 < C 2 , r > p^ 4 and satisfying 



( , 43) M + M + m+ ±^ + ^ + 0)<c Wli , 



with constant C(N) uniform in n and f] 2 ^ h, provided po is sufficiently large. We put 

3/4 

u :— fj, — 11(772), e := ^7—, and 5 := h(r/ 2 ) — h(i(r/ 2 )). Then (I7.42p is equivalent to 

(7.44) P(z/, 5, e; r/ 2 ) := z/ 2 + - e 3/4 ) - + e 2 = 0. 

Pro Pro P« 

Here, we temporarily consider S as independent variable, and a±, a 2 , and /3 are considered 
as functions of v and e, depending on ?7 2 and p n as parameters. It follows from corollary 
12.31 that there exists a neighborhood to of (u, S, e) = (0, 0, 0) such that F = in u if and 
only if 

(7.45) v 2 - 2X 1 (5, e- V2 )u + X 2 {6, e; 772) = 0. 

Here, A 1? X 2 are analytic in 5, e in w and Xi(Q, 0; 772) = A^O, 0; 772) = (the difference 
between ( 17.451) and (I7.44p is that functions X±, X 2 do not depend on v). Moreover, 
it follows from corollary 12.41 and uniformness of our estimates that u can be chosen to 
depend on V only; we also can achieve that uj contains the set |e| < eo, \8\ < 5q, where, eo 
and Sq do not depend on n and r/ 2 G I\ (they depend only on V). We also have uniform 
upper bounds for X%, X 2 and its derivatives. Indeed, we have uniform upper bound 
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( 12.51) . Then analyticity of W (which is equal to v 2 — 2X\(5, e; 772)^ + X 2 (5, e; 772) in our 
case), f)2.5p . and Cauchy's integral formula imply upper bounds for the coefficients X 1; 
X 2 . Now we can solve quadratic equation (17.45p and obtain = X± ± \/X\ — X 2 , or 
^1,2 = h(r/ 2 ) +X X ± y/Xf - X 2 . Thus, recalling that H(£) = r 2 + Z(r) = r 2 + i?" 2 Z(r), 
we deduce that the points 

r 2 + R~ 2 h(V2) + ^" 2 (Xx ± yJxl-X 2 ) = r 2 + h( m ) + R~ 2 {X X ± yj X* - X 2 ) 

are eigenvalues of H(£); notice that, since H(£) is a self-adjoint operator, this implies 
that X\ — X2 > when all the variables take real values. Now, the definition of the 
mapping g implies that there exist two points, a + ,a G Y(£) such that <7(cr±) = 
r 2 + h(r] 2 ) + R~ 2 {Xi ± y/Xf - X 2 ). Lemma EE] implies that a ± G *o(^)- 

Now we recall that S in fact is not an independent parameter, but 5 = h(i] 2 ) — \i(l(t] 2 )) . 
Then our functions X\, X 2 will be analytic in e for |e| < eo, provided ^(772) — h(t(^ 2 ))| < 
5q. Thus, we have proved the following statement: 

Lemma 7.7. There exist positive numbers po, 5q and eo and two functions X\, X 2 , 
Xj = Xj(5,e;r] 2 ), such that Xj are analytic in 5 and e when \S\ < 5o and \e\ < eo and 
the following property is satisfied. Suppose, £ = (£1,^2) £ 2 5 (0) and rj G ^o(^) with 
r/ 2 G h(s), s < min(5 , \C 2 ). Then, there exist two points ct± = ac±(r) 2 ) G ^o(^) such 
that 

(7.46) g(a ± ) = r 2 + h(r, 2 ) + R~ 2 (X X ± y / X 2 -X 2 ), 

3/4 

where Xj = Xj(h(r] 2 ) — h(i(r] 2 )), ^-\r\ 2 ). Moreover, each point v G ^o(^) H l\ can be 
expressed as v = ot±(fi 2 ) for some /j, G T (^). 

Proof. The last statement is the only one we have not proved so far. However, it follows 
from the standard pigeonhole arguments based on the fact that the number of pairs 
(ck+, ck_) is the same as the number of pairs (r] 2 , L(r] 2 )) with r\ 2 G Ii{s). □ 

Remark 7.8. Instead of assuming that n (and thus p n ) is sufficiently large, we prefer to 
assume that po is large enough and prove that estimates hold uniformly for all n (recall 
that p n = 2 n p ). Also, since we always have 5 = h(r] 2 ) — h(i(ry 2 )) depending only on rj 2 , 
we will often skip mentioning the dependence of the functions Xj on the first variable 
and write Xj = Xj(e;r] 2 ). 

Note that Xj and their derivatives enjoy uniform upper bounds and we also have 
Xj(0, 0; 7/ 2 ) = 0. Therefore, by decreasing the values of 5 and e if necessary, we can 
achieve that \X± ± a/ Xf — X 2 \ < C 2 /10 when \5\ < 5 and |e| < e uniformly over rj 2 . 
Now we can fix the value of the parameter s which we used to define the sets h(s) and 
Iq(s): we put s := min(|(5 , \C 2 ). 

Suppose now that 77 G Iq(s). We have: 
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Lemma 7.9. There exist positive numbers po and eo and a function Y = Y(e;r] 2 ), such 
that Y is analytic in e when |e| < e and the following property is satisfied. Suppose, 
£ — (£1162) £ 2 5 (0) and 77 G ^o(£) with r\ 2 G Io(s). Then, there exists a point (3 = 
Pim) e *o(£) such that g(/3) = r 2 + h( V2 ) + R- 2 Y(^;r] 2 ). 

The proof is similar to the above and is even simpler; in fact, this proof essentially is 
equivalent to the proof of lemma 6.1 from [8]. That is why we just give the sketch of the 
proof and make some remarks on uniformness. We start with the formula (17.371) . Now, 
Po is a projector onto one- dimensional subspace corresponding to eigenvalue h(i] 2 ) of B. 
We consider /x such that (cf. (17.391) ) 

\fi-h( V2 )\ <s/3. 

Then operator P(B — /i)P + P is invertible. If necessary we increase po to ensure that 
operator is invertible. Next, we repeat all further arguments from the proof of lemma 
17.71 which are simpler in this case since G M is a scalar- valued function now. We obtain 
that /i is an eigenvalue of Z(r) if and only if 

F := 11(772) - n + -a = 0, 
r 

where a is analytic in \fi — h(r/ 2 )| < s/3, r > p^ 4 and depends on r] 2 as a parameter. It 

also satisfies estimate similar to (I7.43P uniformly in n and r] 2 G Iq. Applying Corollary l2.3l 

(alternatively, we can just use the implicit function theorem) we obtain that F = in 

3/4 3/4 
some neighborhood to of (/i, ^— ) = (h(r] 2 ),0) if and only if /i = h(r/ 2 ) + Y{^-;rj 2 ) for 

3/4 

some analytic function Y . The neighbourhood to contains the set \^—\ = |e| < eo, where 
eo does not depend on n and r\ 2 G Iq. 

As above, we have Y (0; 772) = and Y and its derivatives being uniformly bounded 
(actually, the bound depends on 5 only). Thus, by decreasing e again if necessary, we 
can achieve that \Y\ < s/3 whenever |e| < eo- 

Lemma 7.10. (A) Suppose, i] 2 G h{s). Then we either have oti = 77 and ct 2 = t(r)), or 
ol 2 = Tj and ol\ = l(t]). 

(B) Suppose, r/ 2 G Io(s). Then we have f3 = rj. 

Proof. Suppose, 77 and v are two different points from ^o(^) and v 7^ ^{rf). Suppose for 
definiteness that r(r]) < t(u), i.e. that h(r] 2 ) < h(v 2 ). Then: 

(a) if both 772 and v 2 belong to I , we have h(u 2 ) — h(r/ 2 ) > s, so 

(7.47) h(77 2 ) + y(«; 772) < h(u 2 ) + ¥{■; v 2 ), 

and thus g(fl(r])) < g(f3(v)). 

(b) if both 772 and v 2 belong to I\, we have h(z/ 2 ) — h(7? 2 ) > 3C 2 , so 
(7.48) 

h(77 2 ) + X 1 {- ] 77 2 ) ± V /X 1 2 (.;77 2 )-X 2 (-;772) < h(u 2 ) + X^-) U 2 ) ± yj Xf (• | Z/ 2 ) - X 2 (• \ V 2 ) , 
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and thus g(oc(r))) < g(cc(v)). 

(c) finally, if we have say r\ 2 G h(s) and v 2 G Iq(s), we have h(z/ 2 ) — h(7] 2 ) > 3C 2 , so 



(7.49) h( V2 ) +X 1 (-; V2 ) ± y/X&- 772) - X 2 (-; Tfe) < h(i^) + ¥(■; u 2 ), 

and thus g(ct(rj)) < g((3(u)). 

In all these cases, we have t(rf) < t(u). Now the proof follows from the pigeonhole 
argument. □ 

Corollary 7.11. Let £ = (£ 1; £ 2 ) e -5(0). 

faj Suppose, 77 G ^o(£) wii/i ^2 G 27ien, we either have 



x\- 


-x 2 ) 


/xi 


-x 2 . 


xl- 


-x 2 ) 




-x 2 . 



(7.50) g(rj) = r 2 + % 2 ) + R?{X X + \/ X 2 - X 2 
and 

(7.51) g(L(r 1 ))=r 2 + h(r ]2 ) + R- 2 (X 1 
or 

(7.52) ^ (r7)=r 2 + /l(r?2)+i? -2 (Xi _ 
and 

(7.53) ^(r 7 )) = r 2 + % 2 ) + ^ 2 (X 1 

Suppose, 77 G *o(£) wi£/i ^2 G Io(s). T/ien, we /lave 

3/4 

(7.54) 2(77) = r 2 + % 2 ) + R- 2 Y(^; r ]2 ). 

r 

Remark 7.12. Suppose, 77 and 77' are two different points with the same second coor- 
dinate rj2 G I\. Then we either have both £(77) > £(£(77)) and t(rj') > t(t(r) f )), or both 
£(77) < £(4(77)) and £(77') < t(t(rf)). This shows that we either have fl7T5Til) -f l7T511 or 
fl7.52l) - fl7.53p simultaneously for both 77 and 77'. 

Since the derivative of Y is bounded in {|e| < e }, expression (17.541) is increasing 
function of r (assuming, as we always do, that po is sufficiently large). We denote by 
q = q(r) 2 ) (= q(&)) the value of r which makes the RHS of the equation (I7.54p equal to p 2 . 
Unfortunately, the same argument will not work with expressions (17.501) or (I7.52p (when 
we differentiate the RHS of these formulas, we obtain square root in the denominator). 
It turns out, however, that if we fix r] 2 , the equations 



(7.55) r 2 + hfo) + R~ 2 (X 1 + ^X 2 - X 2 ) = p 2 
and 

(7.56) r 2 + ^(773) + R- 2 (X, - J X 2 - X 2 ) = p 2 
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have exactly one solution each. Indeed, the intermediate value theorem implies that 
there is at least one solution to each equation, and later on in remark 17.131 we will see 
that the total number of solutions of (I7.55P and (I7.56P is at most two. We denote by 
q = q(i] 2 ) (= the value of r which makes the RHS of the relevant equation (I7.50p 

or (17321 equal to p 2 . 

Then similarly to our proof of lemma 14.101 (more precisely, of equation (14.201) ) , we 
obtain the following formula: 

(7.57) voi(A+ n 4 n i) - voi(i- n 4i) = f - \jp 2 - v!)dV2- 

J —a 

Thus, in order to compute vol(A + (lA^^) — vol(v4~ DA^), we need to compute (7(772). 
We will consider the case where 772 G I\ (another case is simpler and can be dealt with 
in the same way). We also assume for defmiteness that formulas (I7.50I) - (I7.51I) are the 
valid ones, so we need to solve equation 



(7.58) r 2 + h( V2 ) + R^iX, + yJXl - X 2 ) = X. 

Thus, (7(772) is the (only) value of r which makes ( 17. 58ft valid, and 5(^(772)) is the (only) 
value of r which solves the following equation: 



(7.59) r 2 + % 2 ) + R- 2 (X 1 - ^ X 2 - X 2 ) = A. 
Now we introduce a new unknown variable 

(7.60) a := - - 1, 

P 

so that 

(7.61) r = p(l + a); 
we also put 

(7.62) ~e:=pT/P- 

Then direct calculations show that (17.581) is equivalent to 

(a 2 + 2a) + p; 3/2 e 2 (/7(77 2 ) + R- 2 X x {e{l + a)' 1 ; 772)) 

(7.63) 



- Pn 3/2 K 2 e 2 ^X 2 Ce(l + a)-\ ^ - X 2 (e(l + o)~\ m ). 
Taking square of the last equality we obtain 

(7.64) W(a, e; 772) := (a + 2) 2 a 2 + i 2 p~ 3/2 w(^ 772) = 0, 

where w is a certain function; the properties of w follow from the properties of Xj. In 
particular, w is analytic in \a\ < 1/2 and |e| < eo/2. Moreover, the bounds for w and its 
derivatives are uniform in n and 772 G 1%. We see that 

iy(0,0;772) = C(0,0;77 2 ) = and W^(0, 0; 772) = 8 ^ 0. 
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Applying again theorem 12. 1\ we obtain that (in the neighborhood of (a, e) = (0,0)) 
W(ct, e; 772) = if and only if 

(7.65) a 2 - 2X 3 (e; 77 2 )a + X 4 (e; r/ 2 ) = 0, 
where Xj(e; 772) are analytic in e and Xj(0; 772) = for j = 3, 4. 

Remark 7.13. Since the equation (17. 65ft has two a-solutions, this implies that the total 
number of solutions of (17.551) and (17.561) is at most two. 

The solutions of ( 17.651) are ai(e,r] 2 ) : = X 3 + a/X| — X 4 and cr 2 (e; 7/2) : = X 3 — 
a/X| — X 4 . Thus, we either have qfa) = p(l + and q(t(j]2)) = p(l + 0^), or 
9(^2) = p(l + ^2) and q{i{r]2)) = p(l + u\); for the sake of definiteness we assume 
the former possibility. 

Put 

T(77 2 ,p) = T„(t7 2 ,p) := q(rj 2 ) + = p(2 + a x + a 2 ) = 2p(l + X 3 ). 

According to Corollary l2.4[ X 3 and consequently p~ l T n is analytic in e for |e| < c(V) with 
some constant c(V) > uniform in n and 772 G ii. Function p _1 T n and its derivatives 
are bounded uniformly in n and 772 G I\. 

Thus, for 772 G h we have p _1 T n is analytic (while 01, a 2 are only algebraic). 

Now assume that 772 G Iq. As usual, this case is similar to the case 772 G Ii, but simpler. 
We solve equation 

3/4 

(7.66) r 2 + h( m ) + R- 2 Y{ — ^ m) = A, 

r 

where Y is analytic in 1/r. Using arguments similar to the first case, we obtain that 

(17.661) has a unique solution r =: 9(772). We define T(r]2,p) = T n (rj2, p) := 2g(?7 2 ) for 

772 G Jo- 
Next, notice that |e| < c(V) is satisfied if r > p^ /4 2(c(y))~ 1 = pT°{p 1 n /2 °c{V)/2)- 1 . 

Thus if we assume that p is large enough, we ensure that the set {r > pt/ 5 } is included 

into domains of analyticity of all our analytic functions. 
According to (I7.57P we have 

voi(i+n4"i)-voi(i-n4t) = 

[ a (q(v2) - \/p 2 - vi)dv2 = \ T(T(p 2 ) - 2 V V-% 2 )^ 2 . 

Note that although 9(772) is not an analytic function of p (it involves square root of 
analytic functions), the function T(r) 2 ) is analytic. 

The proof is almost finished, since the RHS of (I7.67P is analytic in p for sufficiently 
large p (recall that |r7 2 | C a <C Pr/ 3 )- The only remaining thing is to obtain some 
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estimates for coefficients in the analytic expansion of T n (and thus of (17.67!) ). We have 



(7.68) 



T n = p^tjin,^) — . 



j=0 



It easily follows from (17. 63ft . (17.661) and definition of T n that t = 2 and t\ = 0. Since 
P _1 |^n.| — C uniformly in r] 2 for any p > p n , we obtain 



(7.69) 



Un^ 2 )\<C'pf 



with constant C > uniform in n and t] 2 . Substituting it into (I7.67p . we derive 



(7.70) 



1 00 1 

voi(i + n Ail) - n ^Kz) = ~ 2 p E g i( n b 



J=2 



where ej(n) = / (tj(n, — tj{r] 2 ))dr]2- Here, we denoted by t,- coefficients in analytic 

—a 

expansion 



2 \/p 2 -m = p 



We have 
(7.71) 
Next, 



N <4C"p^' /5+1/3 . 



OO oo 

(7.72) p n £ |e-(n)|4 < Cp^ 3 £ p~^ < Cp- 6 ^ 5 + 4 / 3 < Cp^ 

j=6M P' n i=fiM 



i=6A/ 



1 s 



Thus, using (17.701) and 



with constant C > uniform in n. Put ej(n) 
estimates of the coefficients obtained above, we arrive at 

Lemma 7.14. 

GM 

(7.73) voi(A+ n A%) - voi(i- n = E ^("^ + ^ G 4/U 

J'=l 

^|e,(n)|«p^' /5+6/5 . 

Lemma 13.31 now follows after summation over / from lemma 17.141 lemma 16.101 lemma 
14.81 and corollary 14.41 This finishes the proof of lemma 13.31 and, therefore, of theorem 
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